Nucleon Compton scattering in the Dyson-Schwinger approach 
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We analyze the nucleon's Compton scattering amplitude in the Dyson-Schwinger/Faddeev ap- 
proach. We calculate a subset of diagrams that implements the nonperturbative handbag contribu- 
tion as well as all channel resonances. At the quark level, these ingredients are represented by 
the quark Compton vertex whose analytic properties we study in detail. We derive a general form 
for a fermion two-photon vertex that is consistent with its Ward-Takahashi identities and free of 
kinematic singularities, and we relate its transverse part to the onshell nucleon Compton amplitude. 
We solve an inhomogeneous Bethe-Salpeter equation for the quark Compton vertex in rainbow- 
ladder truncation and implement it in the nucleon Compton scattering amplitude. The remaining 
ingredients are the dressed quark propagator and the nucleon's bound-state amplitude which are 
consistently solved from Dyson-Schwinger and covariant Faddeev equations. We verify numerically 
that the resulting quark Compton vertex and nucleon Compton amplitude both reproduce the nyy 
transition form factor when the pion pole in the channel is approached. 
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I. INTRODUCTION 

For several decades, Compton scattering has been an 
important tool to analyze the internal structure of the 
nucleon. At low energies, Compton scattering allows to 
probe the nucleon's ability to change its size and shape 
under the influence of an external electromagnetic field. 
These properties are encoded in the nucleon's polariz- 
abilities for the case of real Compton scattering (RCS) 
and their generalized version for the case of a virtual in- 
coming photon (VCS). Low-energy Compton scattering 
thus provides information on global as well as spatially 
resolved electromagnetic properties of the nucleon and, 
with increasing energy, allows to study the effects of its 
excitation spectrum. In contrast, Compton scattering at 
large energies probes the fundamental quark and gluon 
degrees of freedom. Here, real and virtual photons pro- 
vide important information on the quark and gluon dis- 
tributions and allow to extract the spatial distribution 
of the partons inside the nucleon. Consequently, Comp- 
ton scattering is pursued at many experimental facilities 
around the world including MAMI (Mainz), JLab and 
HI7S at Duke, see e.g. [1-8] for overviews. 

In addition, much can be learned from Compton- 
like processes where two virtual photons are exchanged 
between a charged source and the nucleon [9]. Such 
processes are particularly important in lepton-hadron 
scattering where they contribute to form-factor mea- 
surements and serve to explain the difference between 
polarization-transfer measurements and results using a 
Rosenbluth separation of the proton's electric to mag- 
netic form-factor ratio. Furthermore, the Compton scat- 
tering amplitude in the forward limit contributes to the 
electromagnetic mass shifts in the nucleon system [10] 
and encodes the nucleon structure functions. Its rele- 
vance for the proton radius puzzle has been under recent 
debate as well [11-13]. 



From a theoretical perspective, Compton scattering 
has been studied extensively in effective theories [8, 14- 
18] and from dispersion relations [2, 7], with many im- 
portant results in both approaches. Effective theories are 
a natural tool in the low-energy region, where constraints 
due to chiral symmetry and its breaking pattern play a 
key role. Dispersion relations offer an alternative, reliable 
approach for the analysis of experimental data without 
a-priori limitations for the respective scales. 

While approaches built on explicit quark and gluon 
degrees of freedom are mandatory in the high-energy re- 
gion, their use in the low-energy regime has been limited 
so far. First lattice calculations of polarizabilities have 
been performed, see e.g. [19, 20] and references therein, 
but they still suffer from large pion masses and the omis- 
sion of disconnected diagrams. Model calculations for the 
nucleon polarizabilities are available (see e.g. [21-23] and 
references therein) , but often lack a transparent relation 
to the underlying quantum field theory. 

In this work we begin to explore another approach 
to low-energy Compton scattering. It is based on the 
Dyson-Schwinger equations (DSEs) of motion for QCD's 
Green functions at the quark-gluon level. In combina- 
tion with Bethe-Salpeter (BSEs) and Faddeev equations, 
they provide a comprehensive framework for the study 
of non-perturbative properties of hadrons from QCD at 
all momenta and quark masses [24-27]. In this frame- 
work the masses of baryons have been determined from 
the covariant three-body Faddeev equation [28, 29], and 
electromagnetic as well as axial form factors have been 
calculated [30-32]; see [33] for a brief overview. In a re- 
cent work, we generalized the framework to incorporate 
the coupling of two external currents to the nucleon [34] , 
which opens the possibility to describe a variety of pro- 
cesses such as pion photo- and electroproduction, pion- 
nucleon scattering or Compton scattering in terms of 
their quark and gluon substructure. Here, the impor- 




FIG. 1: Graphical representation of the nucleon Compton scattering amplitude in rainbow-ladder truncation, Eq. (2) with 
Eq. (6). All propagators are dressed, and permutations in the quark lines and symmetrizations of the photon legs are not 
displayed. Diagram (a) corresponds to the term with A^'^G A'^^ and diagrams (b)-(c) to A^" . In terms of Eq. (8), diagram 
(a) illustrates the sum of the first and third lines, diagram (b) the second and diagram (c) the fourth line. The kernels that 
connect the spectator quarks, i.e., the factors (1 — Ka) in Eq. (8), have been absorbed in the hatched amplitudes. Diagram (b) 
contains the IPI quark two-photon vertex and diagram (a) provides the Born parts, i.e., the pure handbag diagrams, to the 
full quark Compton vertex. Note that the electromagnetic current of Eq. (1) in rainbow-ladder truncation has the same form 
as diagram (b) if the two-photon vertex is replaced by the quark-photon vertex. 



tant constraint of electromagnetic gauge invariance dic- 
tates the appearance of well-defined classes of additional 
diagrams besides the Born terms, which naturally ac- 
commodates effects of resonances in s— , u— as well as 
^-channel exchange diagrams. We now apply the general 
framework of Ref. [34] to the specific case of Compton 
scattering. We work out the technical details of the rep- 
resentation of the Compton vertex on the quark level and 
consider first applications in the low-energy region. 

The paper is organized as follows. In Sec. II we 
briefly recapitulate our derivation of the nucleon Comp- 
ton scattering amplitude and discuss the relevant dia- 
grams that appear in the description at the quark-gluon 
level. We isolate the handbag contribution and show that 
it is closely related to the quark Compton vertex. In 
Sees. Ill, IV and V we investigate the general properties 
of a fermion Compton vertex: in Sec. Ill we introduce 
our notation and discuss the Compton-scattering phase 
space; in Sec. IV we present the orthogonal tensor ba- 
sis that we will use in our calculations; and in Sec. V 
we derive a construction for the fermion Compton ver- 
tex that respects its Ward-Takahashi identity (WTI) and 
is free of kinematic singularities. In Sec. VI we present 
first results for the nucleon Compton amplitude and its 
7r'^77 pole contribution, and we conclude in Sec. VII. The 
detailed techniques for solving the inhomogeneous BSE 
for the quark Compton vertex and for computing the 
nucleon Compton amplitude are given in Appendices B 
and C, respectively. 



the Compton scattering amplitude of a baryon: 
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We use a condensed notation where all Dirac, color and 
flavor indices are suppressed and loop-momentum inte- 
grations are implicit. The Lorentz indices and v belong 
to the external currents, and the curly brackets denote 
a symmetrization in these indices. and * / are the 
incoming and outgoing bound-state amplitudes, where a 
bar denotes charge conjugation. Go = S ® S ® S is the 
disconnected product of three dressed quark propagators 
5, and G = Go + Go T Go is the full three-quark (i.e., 
six-point) Green function, with T being the three-quark 
scattering matrix. The quantities A'' and A'"^ are given 
by 
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II. NUCLEON COMPTON SCATTERING 
AMPLITUDE 

A. Construction of the scattering amplitude 

In Ref. [34] we have derived the following combined 
result for the electromagnetic current matrix element and 



Here, and F^"^ are the dressed quark-photon and 
quark two-photon vertices, i^(2) and i^(3) are the two- 
and three-quark irreducible kernels, and K^.y K'^^ are 
those kernels with one or two external photon lines at- 
tached. The subscript 'perm' indicates that each bracket 
has in total three permutations with respect to the quark 
lines. The Compton scattering amplitude is visualized in 
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Fig. 1 for the special case of a rainbow-ladder truncation 
which is given in Eq. (6) below. 

We note that Eqs. (1-2) are completely general. Their 
derivation is based on the formalism of Refs. [35-37] 
which is neither restricted to baryons nor to electromag- 
netic currents and Compton scattering but holds for any 
type of current or hadron. In principle it can also be 
applied to pion electroproduction, nucleon-pion scatter- 
ing, TTTT scattering, etc. If the currents describe pho- 
tons, Eqs. (1-2) are consistent with the hadron's Bethe- 
Salpeter or Faddeev equation so that electromagnetic 
gauge invariance is satisfied. This is only valid as long 
as all underlying ingredients (the quark propagator, the 
two- and three-body kernels and the one- and two-photon 
vertices) are also consistent with each other, in the sense 
that they satisfy Dyson-Schwinger and Bethe-Salpeter 
equations and respect the appropriate Ward-Takahashi 
identities. 

We can extract several features of the scattering am- 
plitude from analyzing the structure of Eq. (2) as fol- 
lows: (a) The scattering amplitude is invariant under 
s — u crossing since the Lorentz indices are symmetrized, 
(b) The first term in (2) with A^'^GA''} generates the 
characteristic handbag diagrams via the three-quark dis- 
connected term in G = Gq + GqT Gq. (c) The same 
term yields also cat's-ears diagram, where the photons 
couple to different quark lines; such diagrams are also 
generated by A'"^ via the second line in Eq. (4). (d) 
The three-quark T— matrix in the interacting term of G 
contains all baryon poles and therefore reproduces all s— 
and u— channel nucleon resonances, starting with the nu- 
cleon's Born term. The scattering amplitude at the re- 
spective pole becomes: 
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(e) Finally, channel meson exchange is implemented as 
well via the quark two-photon vertex F'^'' that appears 
in the first line of Eq. (4) ; it involves a quark-antiquark 
scattering matrix in the i— channel. We will elaborate 
this in the following subsection. 

For our practical calculations we will study Eq. (2) by 
neglecting three-body irreducible interactions (if(3) = 0) 
and using a rainbow-ladder kernel K(^2} — ^rl- That 
kernel describes a dressed gluon exchange with a quark- 
gluon vertex ansatz that is proportional to 7^ and only 
depends on the gluon momentum. Its main advantage 
despite its simplicity is that it satisfies the vector and ax- 
ialvector Ward-Takahashi identities, which in turn guar- 
antee electromagnetic current conservation for form fac- 
tors and the Goldstone nature of the pion in the chi- 
ral limit. Furthermore, it has been shown to describe 
many properties of ground-state pseudoscalar and vec- 
tor mesons as well as nucleon and A baryons reasonably 
well [26-33, 38-42]. Nevertheless, the approximation is 
certainly too simplistic to capture all relevant details of 
the quark-gluon interaction. From a phenomenological 



point of view, the main missing structure comes from 
the pion cloud. Pion-cloud effects are known to be im- 
portant for electromagnetic properties of hadrons such 
as form factors at small momenta or the magnitude of 
polarizabilities. First efforts to include these effects have 
been made in the meson sector [43, 44] along with the 
discussion of other nonperturbative, gluonic corrections 
beyond rainbow-ladder [45-48]. In the baryon sector, 
such extensions have been computationally too demand- 
ing up to now. Consequently, in this exploratory work 
we also restrict ourselves to a rainbow-ladder kernel. 
In rainbow-ladder, the structure of Eqs. (3-4) becomes 
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To maintain all features described above, including 
electromagnetic gauge invariance, we must obtain the 
dressed quark propagator S from its DSE, the quark- 
photon and two-photon vertices from their inhomoge- 
neous BSEs, the nucleon's Faddeev amplitude 4" from 
its covariant Faddeev equation, and the three-body scat- 
tering matrix T from its Dyson equation; all within a 
rainbow-ladder truncation. 

Especially the last part, i.e., calculating the three-body 
scattering matrix T, is a huge endeavor that is not easily 
accomplished with current computational resources. In 
this paper we will focus exclusively on the handbag and 
i— channel structure, i.e., we will ignore the part involv- 
ing the T— matrix as well as all cat's-ears-type diagrams. 
We are aware that electromagnetic gauge invariance can- 
not be realized with this subset of diagrams but only by 
taking into account the full structure of Eqs. (2) and (6). 
Consequently, information on the low-energy behavior of 
the nucleon's scattering amplitude and the generalized 
polarizabilities it contains will be limited, at least in these 
initial studies. 



B. Handbag and channel diagrams 

Our goal in the following is to isolate the nonperturba- 
tive 'handbag' content in the scattering amplitude (2). It 
comprises all diagrams where the two photons couple to 
the same quark line which is, in addition, disconnected 
from the two remaining quarks. For the derivation we 
temporarily suppress all occurrences of the quark propa- 
gator in our notation. In turn, we include the quark-line 
permutations a = 1,2,3 explicitly, so that the rainbow- 
ladder truncated current is given by 



(7) 






FIG. 2: Upper panel: Separation of the fermion Compton ver- 
tex into Born terms and a IPI part. Lower panel: inhomoge- 
neous BSE for the quark Compton vertex in rainbow-ladder 
truncation, where the qq kernel K additionally simplifies to a 
dressed gluon exchange. 



and the combination of Eqs. (2) and (6) becomes 
= */ - ^^«) ri^ (1 - Ka)] 

a 
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r;^ and T'^'^ are the one- and two-photon vertices that 
act on quark a, and Ka is the kernel that connects the 
remaining two quarks. The handbag structure we are 
interested in can only emerge from the first two lines of 
Eq. (8) whereas cat's-ears contributions (a ^ a') appear 
in the third and fourth lines. The three-quark Green 
function becomes in this notation: G = 1 + T, and the 
T— matrix satisfies Dyson's equation: 

3 

T = K{l + T), K = Y,Ka. (9) 

a=l 

To isolate the handbag structure, we extract the two- 
quark contribution with respect to the spectator quark a 
from the three-quark T— matrix: 



T = Tn.+T' . 



(10) 



where T' is the remainder that connects quark a with 
either of the remaining two quarks, or all three quarks 
via iteration of Eq. (9). The last equation entails 

{l-Ka){l+Ta) = l, (11) 

and inserting this together with G = 1 + Ta + T' into 
Eq. (8) yields 



(12) 



plus further diagrams where photons couple to different 
quark lines, including those with T— matrix insertions 
that connect quark a with its companions. Therefore, 
Eq. (12) constitutes the handbag part of the Compton 
scattering amplitude. Comparison with the electromag- 
netic current (7) shows that both equations have an iden- 
tical structure: both include only one factor (l — Ka) and 
the vertices in both equations act on quark a only. The 
quark-photon vertex F'^ in the first expression is replaced 
by a quark Compton vertex 
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in the second, cf. Fig. 2. It is the sum of the one-particle- 
irreducible (IPI) part F^"^ plus a Born contribution Fg"^ 
at the quark level. We will discuss this separation and its 
implications in more detail in Section V. We note that the 
nucleon Compton amplitude from Eq. (12) is diagram- 
matically represented by Fig. 1(b), with the IPI vertex 
replaced by the full quark Compton vertex from Fig. 2. 
The respective diagrams including all kinematic depen- 
dencies are also illustrated in Fig. 10 in the appendix. 

At this point the question remains how the quark 
Compton vertex can be obtained in practice. A suitable 
starting point is the inhomogeneous BSE for the quark- 
photon vertex F^: 



pM 
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where Fq 



Z2 47^ is the tree-level vertex. Go = 5 ® S" 
is now the disconnected product of a dressed quark and 
antiquark propagator, K is the quark-antiquark kernel, 
and the second term involves a momentum-loop integral. 
'Gauging' the quark-photon vertex with an additional 
photon leg yields the IPI part of the quark Compton 
vertex: 
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The second term yields the Born contribution because 

(16) 
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The remainder F^" involves K'^ and does not contribute 
in rainbow-ladder but we keep it for completeness. We 
have therefore 



KGo (Fg"^ +F^''). 



(17) 



Adding the Born term on both sides of the equation and 
using Eq. (13), with T^" := F^' 
inhomogeneous BSE for the full vertex: 



Fp*^, one finds an 
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In rainbow-ladder, the driving term Tq'^ simplifies to 
the Born contribution. The inhomogeneous BSE in that 
case is illustrated in Fig. 2. Eqs. (14) and (18) deter- 
mine the one- and two-photon vertices for a given quark- 
antiquark kernel and quark propagator so that their 
respective Ward-Takahashi identities are automatically 
satisfied, cf. Sec. V. Implementing the vertices in Eqs. (7) 
and (12) yields the nucleon's electromagnetic current and 
the nonperturbative handbag part of the Compton scat- 
tering amplitude in a rainbow-ladder truncation. 

By virtue of the IPI part in the quark Compton vertex, 
Eq. (12) contains more than just the usual quark Born di- 
agrams with two spectator quarks. Using Dyson's equa- 
tion in the quark-antiquark channel, T — K + K GqT, 
we can express the Compton vertex as the projection of 
the full Green function onto the driving term Eq": 



^ 



^-1 ^-pfii' 
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The qq Green function G (or, equivalently, the scattering 
matrix T) includes all meson bound-state poles. If the 
respective meson Bethe-Salpeter wave functions overlap 
with Eq"^, these poles will appear in the quark Compton 
vertex, and consequently also in the Compton scattering 
amplitude of the nucleon. In fact, the quark Compton 
vertex is the sole origin of channel poles in the nu- 
cleon Compton amplitude. The remaining cat's-ears di- 
agrams merely contain p— meson resonances in and 
Q'^ from each quark-photon vertex (which is the projec- 
tion of G onto the tree-level vertex Eg , cf. Eq. (14)), and 
the three-quark scattering matrix produces additional s— 
and M— channel baryon poles. Consequently, the nonper- 
turbative handbag content of Eq. (12) contains also the 
full i— channel meson resonance structure. 

The explicit form of the inhomogeneous BSE (18) for 
the quark Compton vertex and the handbag part of the 
nucleon Compton amplitude (12) are given in in Appen- 
dices B and C, respectively. For the moment we note 
that the nucleon's Compton amplitude J^'^ and the quark 
Compton vertex E^'' have completely identical features. 
They are both fermion two-photon vertices with a rich 
Dirac-Lorentz and momentum structure. Both can be 
decomposed into a Born part and a IPI remainder and 
are partially determined by a Ward-Takahashi identity. 
Therefore, before proceeding with the explicit numerical 
calculation, we will use the next three sections for study- 
ing the general properties of a fermion Compton vertex. 



III. KINEMATICS AND NOTATION 

A. Kinematic variables 

A fermion Compton vertex depends on three indepen- 
dent four-momenta which can be chosen as the four- 
momentum transfer A, the average photon momentum E, 
and the average fermion momentum p, cf. Fig. 3. They 
are related to the incoming and outgoing fermion and 
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FIG. 3: Kinematics in the fermion Compton vertex. 



photon momenta via 

P = \{Pt +Pf) , 
E = i(Q + g'), 

with the inverse relations 



A = Q-Q' ^pj-p. 



(20) 
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(21) 



Depending on the features we want to study, we will 
alternately express the Compton vertex in terms of 
T^"'{p,ll,I^) or T^"'[p,Q,Q'). The first choice is more 
convenient for investigating the general symmetry prop- 
erties of the vertex and for constructing orthogonal tensor 
bases, cf. Sec. IV; it is also highly advantageous for nu- 
merical calculations. It resembles the kinematics of the 
one-photon vertex and nucleon form factors, where A 
would be the photon momentum and S = 0. The second 
choice is more natural for exploring the Ward-Takahashi 
identity for the vertex, and for constructing bases that 
are transverse with respect to both photon momenta and 
free of kinematic singularities; see Sec. V. From any set 
of three independent momenta one can form six Lorentz 
invariants, for example 

[p\Q^,Q'\p-Q,P-Q' ,Q-Q'] 
or {p\ A2, p-E, p-A, E-A}. 

In the case of the nucleon Compton scattering ampli- 
tude, i.e., when the vertex is taken onshell, we rename the 
fermion momenta p, Pi,Pf — > P, Pi, Pf', otherwise all re- 
lations remain unchanged. Here we have the additional 
constraint P^ ~ Pj = —M^, where M is the nucleon 
mass, which entails P-A^O and = _ A'^/i. It 
is most convenient to work with the variables constructed 
from A, E, and P: 



A^ 

UP 



t. 



AP 



(1 + i), P-A = 0, 



(22) 



where t is the (dimensionless) momentum transfer vari- 
able and the limit t = corresponds to forward scatter- 
ing. The remaining variables involve the average photon 
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momentum S and we denote them by X, Y and Z: 



M2 



=:a = tX, P-Y.T=:Y, S-A=:Z. (23) 



Here, a hat denotes a normahzed four-momentum and 
the subscript T refers to a transverse projection with 
respect to the total momentum A, i.e. 



= (,5^'' - A^A'') Y." . (24) 



The Lorentz invariants t and X correspond to radial and 
y , Z to angular hyperspherical variables. Y and Z satisfy 
simple relations under photon crossing and charge conju- 
gation, cf. Sec. IV A. 0(4) symmetry suggests that, kin- 
dred to other Euclidean bound-state calculations where 
hyperspherical angles can often be treated via the first 
few terms in a Chebyshev expansion, the dependencies 
of the Compton amplitude on the angular variables Y 
and Z might be small. In this work we consider only 
spacelike momenta to avoid difficulties associated with 
complex continuations in the nucleon Faddeev amplitude 
and singularity restrictions in the quark propagator; in 
that case we have 



t,X>Q, r, Ze(-l,l) 



(25) 



This is illustrated by the colored regions in Fig. 4, which 
also makes clear that Eqs. (22-23) describe the 'natural' 
variables for Compton scattering in Euclidean kinemat- 
ics. Specifically, we work in the frame where: 
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Notwithstanding that, all subsequent equations will be 
written in a Lorentz-covariant (or, if applicable, Lorentz- 
invariant) manner. This is also mandatory as we will 
frequently switch reference frames when calculating in- 
variant dressing functions of the various ingredients that 
enter the scattering amplitude. 

When considering the quark Compton vertex instead 
of the nucleon Compton amplitude, we must again lift 
the onshell constraints. We can arrange the six (now 
independent) Lorentz invariants of Eqs. (22-23) in two 
sets 



r A2 S2 ^ ^ 
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This distinction will be practical in the context of the 
quark Compton-vertex BSE in App. B, where the mo- 
mentum variables in the set A appear only as external 



parameters whereas those in uj will be solved for dynam- 
ically. Since we replaced the nucleon momentum by an 
offshell fermion momentum, we have now 
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The definitions for t and cr in Eq. (28) anticipate know- 
ledge of the nucleon's mass; however, that requirement 
can be easily removed. 



B. Phase space in nucleon Compton scattering 

Next, we want to investigate the phase space in nucleon 
Compton scattering and the applicability of our methods 
to the relevant kinematic limits. To this end we first 
relate the onshell variables {t, X, Y, Z} of Eqs. (22-23) 
to the incoming and outgoing photon virtualities 



(30) 



4M2 ' 4M2 

and the Mandelstam variables. The latter are obtained 
from the squares of the total momenta 



P, + Q = Pf + Q', 



P 



P-Y = P^-Q' = Pf- 



(31) 



in the s and u channels, see Fig. 5. In analogy to our 
definition of the momentum transfer t, it is convenient to 
define modified Mandelstam variables s and u that are 
dimensionless and positive in the physical region: 

(^^-(1+.,. (32, 



M2 



s, t and u are related to the conventional Mandelstam 
variables s, t, u by 

i^-iM'^t, S = M2(1 + s), il = M^{1 + u), (33) 

where the extra minus signs in Eqs. (32) and (33) orig- 
inate from the Euclidean description. The limits s = 
and u = correspond to the location of the nucleon pole 
in each channel, and nucleon resonances in the Compton 
amplitude appear at s > or m > 0. Photon cross- 
ing symmetry (s O u) is expressed by the dimensionless 
crossing variable v, conventionally defined as 



S • P 



s — u 
4Af2 



s — u 
4 



(34) 



The two sets of variables {t, t, t', v} and {t, X, Y, Z} are 
then related via 



^, \=^-{l+X±2VXZ 



= t{l- X)±2iy, 



(35) 
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with V = yfXy/Z^ - l^t{l + t)Y, 

2(r + r0 
^=^t ^' 



or vice versa: 



Y = 



Z = 



t 



l+t y/{T + t' - <)2 - 4tt' ' 



(36) 



Y plays the role of the crossing variable since it is pro- 
portional to V, and Z ~ r — r' is the skewness variable. 
As required, the Mandelstam sum s + u = 2t{l — X) leads 
to Q'^ + Q''^ + s + i+u = 2A'P. 

We can now investigate various kinematic limits in 
the phase space defined by the variables {t,T,T',v} or 
{t,X,Y,Z}, cf. Fig. 4: 

• Real Compton scattering (r = r' = 0) implies 
Z = and X — —1. The remaining independent 
variables t and Y constitute the Mandelstam plane, 
with {s, u} = 2t± 2^yt{l + t)Y. The physical re- 
gion corresponds to i > and Y > 1, with t = 
forward scattering and Y — 1 backward scattering. 
The nucleon polarizabilities are defined in the real 
forward limit (t = 0). 

• Virtual Compton scattering (t' — 0): here only 
the outgoing photon is real. The skewness variable 
becomes a function of X: 



Z = 



1 + X 



2VX 



> 1. 



(37) 



and the remaining variables are t, X and Y, with 
r^^il + X), ,.^±^^^t{l + t)Y. (38) 

Timelike Compton scattering (t = 0) is the 

reverse process with an incoming real photon; 

Z^-'-±£<-l, r'^\{l + X). (39) 
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-1 




/ i ' 

1 



wcs 



X=1 

QQ'=o 



FIG. 4: Sketch of the Compton scattering phase space in 
the variables r and r' at fixed t {upper panel) and in the 
variables X and Z for arbitrary t ^ {lower panel). The 
abbreviations RCS, VCS, WCS and TCS denote real, virtual, 
doubly virtual and timelike Compton scattering. The colored 
area describes the spacelike region. 



• The generalized polarizabilities are defined at 
the intersection of the VCS condition Q'^ — with 
the limit Q ■ Q' = 0, which implies Z = X = 1. 
In this case, the spacelike region reduces to = 
since the dependence on Y decouples, cf. Eq. (35), 
and the only remaining variable is t = r. This is 
the Born limit s = u = which exhibits the onshell 
nucleon pole in both s and u channels. 

• Forward WCS (i = 0, r = r'): the WCS limit 
corresponds to the line t = t' <^ Z = in Fig. 4. In 
the forward limit t = the spacelike region shrinks 
to the WCS line since the dependence on Z decou- 
ples. For forward scattering we must work with the 
variable cr = tX from Eq. (23) since X is no longer 



well-defined. Because of t = and Z = the four- 
vector A'^ vanishes, and hence — Q'^ = E^. 
The remaining independent variables are a and Y 
(or, equivalently, t and xb) with 



T = r = 



(7 

4 ' 



{"J 






-A 











where the Bjorken scaling variable is 



Xb 



2P- S 



(7 

2^ 



(40) 



(41) 



In the physical region Y must become imaginary; 
we have chosen sign conventions so that xb — 0"*" 
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corresponds to y — ^ +ioo. The forward VVCS 
amplitude encodes the nucleon structure functions 
/t, /l, 9tt and gLT- 

It is clear from these considerations that many of the 
interesting physical applications with available experi- 
mental data correspond to phase space regions outside 
the 'spacelike' domain that is defined by Eq. (25) and 
illustrated in Fig. 4. Directly accessible are: 

• the generalized polarizabilities at ^' = 0, including 
the usual polarizabilities in the real forward limit, 

• forward VVCS iov < \Q\/M ^ xg > V^, e.g. 
through a Chebyshev expansion to complex Y, 

• and two-photon exchange processes with a charged 
source where the full spacelike region is integrated 
over. 

Access to the remaining limits (real and deeply virtual 
Compton scattering, Bjorken regime) is in principle pos- 
sible but requires higher numerical effort and sophisti- 
cated methods for dealing with complex singularities. For 
the time being, we will restrict ourselves to those appli- 
cations which can be obtained from the spacelike calcu- 
lation and, if necessary, perform extrapolations to the 
inaccessible phase space domains. 



IV. TENSOR BASIS FOR A FERMION 
COMPTON VERTEX 

In this section we investigate the general properties of 
a fermion Compton vertex and construct a simple, or- 
thonormal tensor basis that is useful for numerical cal- 
culations. The analysis can be applied to the nucleon's 
Compton amplitude as well as the quark Compton ver- 
tex. For the former, the incoming and outgoing nucleons 
are onshell which simplifies the structure considerably. 
We begin with the analysis of the general offshell case 
and return to the onshell simplification in Sec. IV B. 



A. General Compton vertex 

We start by constructing the general Lorentz-Dirac 
basis for a fermion Compton vertex r^''(p, E,A). The 
vertex is a matrix in Dirac space and depends on three 
four-vectors p, S and A, cf. Fig. 3 for the momentum 
routing. We work with the kinematic variables defined 
in Eq. (28): the photon variables t, a and Z and the 
fermion variables p^, z and y. 

To facilitate the analysis we want to express the vertex 
in an orthonormal basis. We can achieve this by orthog- 
onalizing the three momenta A, S and p, so that 



d'' := A^ , 



Pt 



(42) 



are three covariant unit vectors. The hats denote nor- 
malized four-vectors, and the subscript T stands for a 



transverse projection with respect to the external mo- 
mentum A, e.g.: 



= = (^^'^ - A^A"^) Yy . 



(43) 



In addition, the subscript t indicates a transverse projec- 
tion with respect to both A and Sj^, i.e.: 



Pt - Pt v2 



T 



(44) 



The resulting vectors r, s and d satisfy (P = = = 1 
and d ■ s ^ d ■ r ^ s ■ r ~ 0. 

Furthermore, instead of using or some transverse 
combination thereof, the dependence of the vertex on 
three independent momenta allows us to work with the 
axialvector 



(45) 



as the remaining independent four-vector, v must ap- 
pear in combination with 75 to transform like a vector 
under parity, which allows to construct a complete and 
orthogonal tensor basis from all possible combinations of 

75^;^ r^ s^ 1, f, ^, 4, (46) 

with ^ = 75 / ^ 4. For example, we can now write 

7^ = w^i^ + r^'f + s^'^ + ■ (47) 

The structure S'^'^ no longer appears in the basis since it 
is linearly dependent: 

S^"' = v^v" r'^T" -f s^s" ^- d^d"" . (48) 

Similarly, the commutator [7^,7'^] can be expressed as 

+ ^ ^ (w^s'' - s^v") + f4 {r^d" - d^r") (49) 
+ ^4 [v^d" - df'v") + ^4 [s^d" - d^'s") , 

and one can write 

In the Lorentz frame defined by Eqs. (26) and (29), u, r, 
s, and d become the Euclidean unit vectors in R"*: 



, d = 



(51) 

and the covariant relations (47), (48) and (50) can be 
read off immediately. 

The vertex r^''(p, S],A) consists of 128 independent 
Lorcntz-Dirac basis elements. With the relations above 
it is straightforward to construct a basis whose elements 
are 
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ST 


BC 








ST 


BC 


xr 

xr 




+ • 

+ o 


++ 
++ 




xr 

-^10 




+ s^r" 
+ s^i;" 


+ o 
+ o 


-+ 

+- 


xr 
xr 


Li 1^ U f 


+ o 

— o 


++ 




-^11 
-^12 




- s^r" 

- s^u" 


— o 

— o 


+- 
-+ 


x^ 

xr 


11, U 1 HI lh> 

s + d^d 


+ • 

+ o 


++ 
++ 




X'l3 


r^d" 
v^d 


+ d^r" 
+ d^u 


+ o 
+ o 


+- 
- + 


xr 
xr 


s^d" + d^s" 
s"d'' - d^s" 


+ o 

— o 


++ 




xr 

-^16 


r^'d-' 
v^'d" 


- d^r" 

- d^u" 


— o 

— o 


-+ 

+- 





BC 


T\ 


1 


^ — \- 


T2 


/ 


++ 


T3 




- + 


T4 


4 


+- 


rs 


ft 




re 


f4 


++ 


T7 




- + 


rs 


fH 


- + 



TABLE I: Tensor basis for a fermion Compton vertex. The labels S and T indicate whether the respective basis element is 
symmetric or antisymmetric in the Lorentz indices and whether it has a Lorentz trace or not. B and C denote the Bose and 
charge-conjugation symmetries for each element. The full basis is constructed according to Eq. (52). All basis elements must 
be additionally equipped with a factor ^ to ensure orthonormality in the sense of Eq. (59). 



• orthonormal, 

• symmetric or antisymmetric in the Lorentz indices, 

• have a non- vanishing Lorentz trace or are traceless, 

• have definite transformation properties under pho- 
ton crossing and charge conjugation, 

• fully factorize the Lorentz and Dirac structure. 

The result is collected in Table I and involves 16 Lorentz 
structures Xr (« = 1 • • • 16) and eight Dirac structures 
Ta (a = 1 ... 8). The general basis element is constructed 
as: 



f::(r,s,d)=xr(r,.,d) 



1 

75 



ra(r, s,d) 



(52) 



where the bracket means that the negative-parity ele- 
ments which involve only one instance of the axialvec- 
tor (i.e., those with i = 3,4,10,12,14,16) must be 
equipped with a factor 75 to restore positive parity. The 
decomposition of the fermion Compton vertex can then 
be written as 



16 



(53) 



=1 0=1 



where the 128 dressing functions fi a depend on the six 
Lorentz- invariants that appear in the sets uj and A, and 
the basis elements depend on the three orthonormal mo- 
menta r, s, d. The advantage of this construction is 
that Lorentz-invariant relations between the basis ele- 
ments become very simple: all inherent scalar products 
are either zero or one, and Lorentz traces can be conve- 
niently separated from Dirac traces. 

The column labels S and T in Table I define orthog- 
onal subspaces with respect to the Lorentz indices: the 



elements with S = ± are symmetric or antisymmetric in 
the Lorentz indices, those with T = • have a nonvanish- 
ing Lorentz trace, and the ones with T = o are traceless. 
We show in App. B 3 that the three possible {S, T} com- 
binations define different subspaces which decouple in the 
inhomogeneous BSE for the quark Compton vertex. In 
fact, all operations on the Lorentz indices and/or exter- 
nal momenta s'^ and leave the relative-momentum in- 
tegration unchanged, and the basis in Table I leads to the 
maximum number of decoupled equations, namely 11. 

The remaining columns labeled B and C correspond 
to the Bose and charge-conjugation symmetry properties 
of the basis elements. The full vertex must be Bose- 
symmetric (i.e., invariant under photon crossing) and in- 
variant under charge conjugation: 

r'^''(p,i],A) = r^rp,-s,A), (54) 

f'^''(p,S,A) = Cr'''^(-p,-S,-A)^C'^ 

= r^''(p,E,-A). (55) 

The basis elements in Table I are either symmetric or 
antisymmetric under either of these operations. (For ex- 
ample: the element is Bose-antisymmetric but sym- 
metric under charge conjugation.) The symmetry of the 
full basis element t/*^ of Eq. (52) is the combination of 
the individual symmetries for Xr and (e.g., r^g is 
both Bose- and charge-conjugation symmetric). An ex- 
ception are the opposite-parity elements which carry an 
additional factor 75 : when combined with any of the spin 
structures a = 2,3,4,8, they switch C— parity because 
75 has to be permuted back through an odd number of 7 
matrices. 

Furthermore, the momentum variables y and Z from 
Eq. (28) flip their sign under photon crossing, and z and 
Z switch sign under charge conjugation. The variables 
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, t and a are invariant under both. The invariance of 
the full Compton vertex under both operations implies 
that each basis element in Table I with symmetry BC 
must be matched by an appropriate dressing function 
fi,a with the same symmetry. This yields the following 
symmetry properties for the dressing functions fi a for 
given B,C e {±1}: 

/(z, y, -Z) = B f{z, -V, Z) = C f{-z, y, Z) , (56) 

where we omitted the labels z, a and the momentum de- 
pendencies in p^, t and a. Eq. (56) entails that the 
dressing functions can be factorized into angular pref- 
actors which carry the symmetry and remainders which 
are completely symmetric in all variables. For each BC 
combination there are two possible choices for these pref- 
actors: 



BC = ++ 4^ 1 or yzZ , 

BC = — h ^ y or zZ , 

BC = H — 4^ z or yZ , 

BC = 4^ Z or yz . 



(57) 



Finally, the basis has the advantage of being orthonor- 
mal. The orthonormality relation depends on the charge- 
conjugated basis elements, defined by 



Tr^T 



':^ir,s,d) = Cri^^:i-r,-s,^dyC 



= Xr(r,s,d)CT„(-r,- 

= Xf''(r, s,d) fa{r,s,d) 
It is given by 



-dye 



1 

75 



(58) 



1 

75 



= 5ij \ Tr {fa n} = 5ij Sab 





■ 1 ■ 




■ 1 ■ 












/•I 






i 







(59) 



Since the Xj are orthogonal by themselves, Xf'Xj'^ — 6ij, 
the product of the square brackets is the unit matrix. 

As the basis vectors are normalized to unity and con- 
tain transverse projectors with respect to A and E^, one 
will encounter kinematic constraints in the limits — 0, 
E|, = and = 0. This is just the boundary of the 
spacelike region where any of the variables t, a or 
vanishes or Z, z or y become ±1. Hence, the dressing 
functions in those limits are interrelated and/or can have 
kinematical zeros. This is not a problem in the practical 
numerical calculation as long as we avoid to work in these 
exact kinematical limits but rather approach them from 
the interior of the spacelike region. Nevertheless, it may 
affect the numerical accuracy in their vicinity. 

So far we have not worked out the transversality of the 
photon at the level of the basis elements of the Comp- 
ton vertex. This is an important issue since electromag- 
netic gauge invariance entails that the onshell nucleon 



Compton amplitude is fully transverse, i.e., transverse 
with respect to Q' in the index fi and Q in the index 
v. Moreover, also the quark Compton vertex can be 
split into a piece that satisfies its Ward-Takahashi iden- 
tity (containing longitudinal and transverse parts) plus 
a remainder that is purely transverse, cf. Sec. V. Since 
the complete dynamics of the system is encoded in the 
transverse pieces, it is in principle sufficient to work with 
a transverse basis only. The distinction between trans- 
verse and longitudinal remains intact when the quark 
Compton vertex is implemented in the nucleon Compton 
diagrams, and electromagnetic gauge invariance requires 
that the sum of all longitudinal diagrams must cancel at 
the nucleon level. 

In Sec. VE we will construct a fully transverse basis 
that is free of kinematic singularities; unfortunately that 
basis is no longer orthogonal and therefore cumbersome 
in the numerical calculation. Here we will pursue a dif- 
ferent goal and construct a transverse basis from Table I. 
Since those basis elements factorize the Dirac from the 
Lorentz structure, it is sufficient to analyze the transver- 
sality relations at the level of the 16 Lorentz structures. 
The transverse part satisfies the conditions Q'^ F'^'^ = 
and F^"^ = which can be worked out to construct 
a fully transverse, orthonormal basis that consists of 9 
elements, see App. D: 

Yi = Xi , 

Y2 = X2 , 

(1 - X)(X5 + Xs) - 2b (Xs - aXy - bXe) 
\/2nin2 

Y4 = X4 , 
Y5 = X3 , 
1 



Y., = 



Ye = — (Xn + aX9-6Xi3) 
Yy = ^ (Xg + aXn~ bX,,) - 
Y8 = ^(Xi2 + aXio-feXi4), 



(60) 



Yg = -L (Xio + aXi2 - 5Xi6) - 4^ 



We have abbreviated 



a = VXZ, b = VX ^/l- Z^ 



(61) 



and 



ni 



1 + X . 



n2 = ni 



4a2 
ni 



The new basis elements Y^ can now have kinematic sin- 
gularities if ni = or 712 = 0. From Eq. (37) we see that 
^2 = in virtual and timelike Compton scattering while 
rii vanishes for real Compton scattering, cf. Fig. 4. The 
respective dressing functions must have nodes in these 
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BC 


Yi 


++ 


Y2 


++ 


Y3 


++ 


Yt^ 


++ 


Y475 


++ 


Ys 75^ 


++ 





BC 




BC 


Yi^ 


- + 




Ye 


+- 


Y2^ 


- + 




Y8 75,^ 


+- 




- + 




Y975 


+- 


Y7 


- + 




Ye^ 




Y8 75 


- + 




Y4 75,^ 




Yg 75 ^ 


-+ 




Y575 





TABLE 11: Transverse, orthonormal tensor basis for the on- 
shell nucleon Compton amplitude, together with the photon- 
crossing and charge-conjugation symmetries for each element. 
The tensor structures Yi correspond to the dressing functions 
Fi in Eq. (67) and the structures Yi ^ to the functions d. 



limits so that the singular basis elements decouple. For 
example, the elements Y3 7 g would decouple in the VCS 
limit, so that only six independent (Lorentz) basis ele- 
ments survive. 

The BC symmetries for the Y^ are inherited from the 
Xj; note that the variable a switches sign both under 
charge conjugation and photon crossing. The full trans- 
verse basis is then constructed in complete analogy to 
Eqs. (52)-(53), namely by multiplying the Y^ with the 
Dirac structures from Table I and with appropriate 
75 insertions for the elements i = 4,5,8,9. Therefore, 
we arrive in total at 9 x 8 = 72 independent transverse 
Lorentz-Dirac elements. For completeness, the remain- 
ing 7 X 8 = 56 longitudinal and transverse-longitudinal 
basis elements Y^g are worked out in App. D. The full 



vertex is expressed as 



16 8 



i—l a—1 

with 75 insertions for i = 4, 5, 8, 9, 15, 16 



1 

75 



Ta{r,s,d) , 
(62) 



B. Nucleon Compton amplitude 

Next, we apply our general considerations to the spe- 
cial case of the nucleon's Compton amplitude. To this 
end we replace the relative momentum p again by the 
average nucleon momentum P and the variable y by Y . 
The onshell constraints Pf ~ Pj = ~AP entail 



P' 



-M^{l + t), z = P-A = 0, (63) 



so that there are four instead of six independent Lorentz- 
invariants. P is now already transverse to A, i.e., Pt — 
P. Sandwiching the vertex between onshell spinors is 
equivalent to applying the positive-energy projectors 



f 



from left and right, so that the matrix- valued nucleon 
Compton amplitude has the structure 



J'"'(P, S, A) = A( ri^'iP, E, A) A*^ 



Eq. (63) 



(65) 



Acting with the projectors on the Dirac basis elements 
Ta(r, s, d) leaves only two independent structures, namely 
1 and ^ (or, equivalently, 1 and ?5), whereas the remain- 
ing ones are either zero or become linearly dependent via 
the Dirac equations 



A^ff = iMA{ , AX = iMA\ 



(66) 



From Table I, this leaves 16 x 2 = 32 independent ten- 
sor structures. In contrast to a general fermion Compton 
vertex, the onshell nucleon Compton amplitude is fully 
transverse, cf. Eq. (105) below. The transversality condi- 
tion reduces this number to 9 x 2 = 18, and it is sufficient 
to work with the Y^ given in Eq. (60) and Table II. The 
resulting Compton amplitude assumes the form 



J''''(P,E,A) ^^Yf'^A^; 



1=1 



1 

75 



(F, + G,^)A;, (67) 



where the Fi and Gi are functions of X, Y and Z . All 
other symmetry relations are maintained, and Eq. (56) 
simplifies to 



f{Y,-Z) = Bf{-Y,Z) = Cf{Y,Z) 



(68) 



BC = ++ 




1, 


BC = -+ 




Y, 


BC = +- 




YZ, 


BC = — 




Z . 



(64) 



for / e {Fi,Gi}. The possible symmetry prefactors of 
Eq. (57) become: 



(69) 



These relations have also been studied in the context of 
dispersion relations to determine the low-energy behav- 
ior of the individual dressing functions of the Compton 
amplitude [2, 49-51]. (Note that from Eqs. (34) and (36): 
Z ^ -Q'^ andY ^ P = P -Q = P ■ Q'). 



V. WARD-TAKAHASHI IDENTITY AND 
TRANSVERSALITY 

The basis construction of the previous section is conve- 
nient for studying the general properties of the Compton 
vertex and also for the numerical implementation. How- 
ever, it is less useful for working out the properties re- 
lated to electromagnetic gauge invariance, i.e., the Ward- 
Takahashi identity, the transversality conditions and re- 
lated analyticity constraints. We will deal with these 
issues in the following and decompose the Compton ver- 
tex into a 'gauge part', i.e., the contribution that is fixed 
by its Ward-Takahashi identity, and further transverse 
pieces that are free of kinematic singularities and con- 
straints. 
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A. Fermion-photon vertex 

We start with a discussion of the fermion-photon ver- 
tex as it provides the template for the two-photon case. 
It satisfies the Ward-Takahashi identity 

Q''r^'{k,Q) = S-\k+)-S-^{k^), (70) 

where Q is the photon momentum, k is the relative mo- 
mentum of the quark, and k± = k ± Q/2 are the quark 
momenta. The inverse dressed quark propagator reads 



S-\k) = 4 A{k^) + B{k^) , 



(71) 



and the renormahzation-point independent mass func- 
tion of the fermion is given by M{k'^) = B{k^)/A{k^). 
Eq. (70) is solved by the Ball-Chiu vertex [52] 



where the functions 

SA(fc,Q) 
AA(fc,Q) := 



(72) 



As(fc,Q) 



A{kl) 


+ A(fc2) 




2 


A{kl) 


-A{kV) 


k\ 




B(kX) 





(73) 



kl-k^_ 



are completely determined by the dressed fermion prop- 
agator and free of kinematic singularities. 

The full vertex is then the sum of the Ball-Chiu part 
and a transverse piece that is not constrained by the 
WTI: 



r^(fc,g) =rg^(fc,Q) 



r^(A:,Q) 



(74) 



ri^ consists of eight independent tensor structures. An- 
alyticity at vanishing photon momentum requires to 
vanish in the limit = 0, either via appropriate mo- 
mentum dependencies of the basis elements, vanishing 
dressing functions, or kinematic relations between the 
dressing functions in that limit. In order to find eight 
kinematically independent dressing functions, we want 
to express in a basis that is free of kinematic singu- 
larities and 'minimal' with respect to its powers in the 
photon momentum. Since the construction of the two- 
photon vertex is closely related to the one-photon case, 
we illustrate the problem here in detail. 

The general fermion-photon vertex with quantum 



numbers J 



PC ^ 1 — 



1 vertex consists of 12 tensor struc- 



tures which can be chosen as 



(+) r 

(-) [rJ] 
(+) 

(+) b'',l(, 



(+) (-) 

(+) k^^ (-) g'^^ 

(-) k^Q (+) Q^g 

(+) k'^[^,m (-) Qn^,m- 



(75) 



To ensure definite charge-conjugation symmetry (indi- 
cated by the signs in the brackets) we have used the 



commutator for the product of two 7 matrices and the 
totally antisymmetric combination 

[A, B, C] [A, B]C+ [B, C]A+ [C, A] B (76) 

for three 7 matrices. If the odd basis tensors are multi- 
plied with a factor k ■ Q, the full vertex satisfies 

r^(fc, Q) = cr^(-fc, -Qfc^ = -Tf'ik, -g) (77) 

with scalar dressing functions that are even in A; • g. 

The transverse part of the vertex consists of eight 
tensor structures that are constructed from Eq. (75). 
The two elements [7'', (^] and [7'^ , ^, (^] are transverse by 
themselves. In principle one could apply the transverse 
projector 



(78) 



to the remaining elements from the first two columns of 
Eq. (75) to obtain the basis decomposition 



-i^T = 5i7t + 92 k-Q f [7^, ^] 

+ kt^ {igb + 9(>^ + 97k-QQ + gal 

where 



11 rTiUl> U 

7t = 7 : 



(79) 



(80) 



We have attached prefactors so that the scalar dressing 
functions gi{k'^ , k ■ Q, Q^) are even in k ■ Q and real for 
fc^ > 0, g^ € R. However, since the projector (78) con- 
tains a kinematic singularity at 0, the resulting 
dressing functions are kinematically dependent: the four 
combinations 

gi + {k- Qfgi , .92 - g& , 35 , ge (81) 

must vanish with for — >■ 0. Instead of the pro- 
jector (78) one could equally apply Tq' which has 
no kinematic singularity; unfortunately this overcompen- 
sates the problem since gi, g2^ gr, ga do not need to vanish 
individually when goes to zero. 

A basis decomposition where all dressing functions are 
truly kinematically independent is given by [53-55] 



-iK^hQ^l'T + f2k 



+ /3i[7^^] + /4H7^^:^] 
+ ih 0" fc^ + h K $ 

^ hk-Q{k-Q-i^ -k^% 
+ hW-Ql^~k^(liM 
It satisfies the requirements of Eq. (81) since 

/iQ' = .91 + (A: • Q) V , /sQ' =35, 

/2 = .92 - 58 , /e = 36 , 

/a = .93 , -/t = 37 , 

/4 = 34 , /s = 38 • 



(82) 



(83) 
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Apart from global factors k-Q, the four tensor structures 
corresponding to /3,4,7,8 are linear and the remaining four 
are quadratic in the photon momentum. 

The question remains whether Eq. (82) can be ob- 
tained from a systematic construction principle. To this 
end we define the quantities 



75 e'"'"'^a"b'^ , 



(84) 



with a'^,bf^ € {fc*", Q^}. They are both regular in the 
limits a — !■ or 6 



0. is transverse to a'^ and b"^, 



■ Mf u-^ — 



(85) 



whereas e(^^ is transverse to a and b in both Lorentz in- 



dices. The usual transverse projectors can thus be writ- 
ten as Tj^^ = t^^/Q^ and T^r = t^Q,Q,/Q''. 

With the help of these definitions one can generate the 
basis (82) as follows. Take the four tensor structures that 
are independent of the photon momentum: 



7^ [rj], k\ k'^^. 



(86) 



Contract them with tgg 



''Qk 



and Sq^ to generate eight 
transverse basis elements that are kinematically indepen- 
dent and linear or quadratic in the four- momentum Q^: 




(87) 



k■Q-f^'-k^'Q 1 
[k-Qr-kf^QJ] j ' 

Instead of using and eg^, , one could contract the four 
elements in Eq. (86) also with i^"^ ^ Sf"" - -y^'Q" and 
use commutators where necessary. However, this does 
not generate any new elements: 

K,k-'J]^-t'i,l[rJ]. 

Finally, attach appropriate factors k-Q to ensure charge- 
conjugation invariance of the dressing functions. 

We will henceforth use Eq. (82) as our reference basis 
for the transverse part of the fermion-photon vertex. We 
write it in a compact way: 



-QQ k-Q l[ 

i[7^^], 



— f'"' ik" 

— Lqq f , 

= C^■Q7^ 



(89) 



.A' _ 1 



~ ^Qfe 2 [7'^J V\ 



The full vertex is thus given by Eq. (74), with the trans- 
verse part 



ir^(fc,g) 



fi{k , k ■ 



')rl^{k,Q). (90) 



The dimensionful dressing functions fi{k^, k ■ Q, Q^) are 
again even in fc • Q. They are kinematically independent 
and can remain constant at vanishing photon momen- 
tum. The basis (89) is essentially identical to Eq. (A. 8) 
in Ref. [53] and Eq. (A2) in Ref. [55]. The relations be- 
tween our rj^ and the transverse tensor structures in 
those papers are 



T2 = -\k-QTi. 



T5 =Ti, 

^k-QTe. 



T6 ^\T2 



T7 - 
T8 ^ hTj . 



(91) 



The dressing functions associated with T3 and T4 con- 
tribute to the onshell anomalous magnetic moment, 
cf. Ref. [48] and Eq. (96) below, and ry constitutes the 
transverse part of the Curtis-Pennington vertex [56] . 

Finally, to obtain a connection with the nucleon's on- 
shell current, we investigate the limit where the incoming 
and outgoing fermion lines are taken on the mass shell, 
i.e. 



kj. — —m or 



m2-QV4, k-Q 
The onshell vertex 

J^(A:,Q) = A^r^(fc,Q)A; 



0. 



Eq. (92) 



(92) 



(93) 



is sandwiched between Dirac spinors that are eigenvec- 
tors of the positive-energy projectors 



Ai 



A+(fc+), 
A+(fc_), 



A+(P) 



1 



(94) 



By virtue of the projectors, only two of the basis elements 
in Eq. (89) remain independent, and the vertex can be 
written in the standard form 

J^{k, Q) = iKi {f^ 7^ + 2 [7^ i^]) AV , (95) 

where Fi , F2 are dimensionless functions of Q^ only. Via 
Eq. (74) they consist of Ball-Chiu parts and transverse 
components which are related to the functions S/i, /S.A^ 
As and fj in the onshell limit: 

Fi(g2) = A{-m^) + 2m [B' {-m^) - mA'{~m?)\ 



fi-m (/s + mfe) 



fi - mfs 



Eq. (92) 



mi 

2m 



2 



2 



Eq. (92) 



(96) 
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FIG. 5: Separation of the fermion Compton vertex into Born terms and a IPI part. 



Here we have exploited the fact that E^j, A^, 
simphfy on the mass shell to A{—m'^), A'{—m?) and 
B'{—m?), respectively. Since m is defined via Eq. (92), 
i.e., m? = —k\ — — fci, all relations hold even in the 
case where the fermion propagator has no timelike pole 
(although in that case Eq. (93) cannot be derived from a 
pole condition but is merely a definition). On the other 
hand, if a mass shell exists where the propagator behaves 
like a free particle, we have in addition M(—w?) = to, 
M'{—m?) — and A{—m?) = 1, and the square brackets 
B' — mA' in Eq. (96) vanish. The dressing functions /a 
and /4 contribute to the fermion's anomalous magnetic 
moment. 



B. Compton vertex, Born terms and WTI 

We proceed by generalizing our findings to the two- 
photon case. The IPI part of the fermion Compton 
vertex satisfies the following Ward-Takahashi identity, 
see [57] and references therein: 



Q"- r^''(p, E, A) = r-(p+, Q) - v^ip-.Q) , 
T^^ip, E, A) = r^(p+, -Q') - T^{p-,-Q') 



(97) 



Here, V^^k, Q) is the dressed fermion-photon vertex that 
depends on a relative momentum k and a total photon 
momentum Q. The relative momenta that appear in the 
fermion-photon vertices are 



p1 



2, .-2' 

and p is the relative momentum in the Compton ver- 
tex. The second equation in Eq. (97) follows from the 
first one via Bose symmetry (replace ^ ^ Q ^ —Q' 
and thus E o -E, and use Eq. (54)). The full Comp- 
ton vertex f F^" + T^'^ is the sum of the IPI part 
and a one-particle reducible, Bose-symmetric and charge- 
conjugation invariant Born piece, cf. Fig. 5: 

rr = - rnp+, -Q') s{p + E) T^p+,Q) 

-F-(p7,g)5(p-E)F^(p,-,-Q'), 

where S{p) is the dressed fermion propagator. 

To obtain the WTI for the full vertex F'"^, we insert 
the WTI for the fermion-photon vertex, 

Q^T^'{k,Q)^S-\k+)-S~\k^) (100) 



(99) 



with k±^k± Q/2, in Eq. (99): 
Q"^F^''(p,E,A) = 

= [S-\p+^)-S-\p + E)] S{p + E) F''(p+, g) 
+ F-(p7, 0) S{p - E) [S-\p - E) - S-\p - f 



We expressed the photon momenta Q and Q' by the av- 
erage photon momentum E and the momentum transfer 
A. Adding Eq. (97) yields the result 



Q'^f'^''(p,E,A) = 

+ S-\p+^)S{p + ^)T''ip+,Q) 
-F''(p7,Q)5(p-E)^-i(p-f). 



(101) 



We can furthermore attach quark propagators from the 
left and right to obtain 



5(p+f)F^''(p,E,A)5(p-f) 
= 5(p+E)F''(p+,g)5(p-|) 
-S{p+^)r''{pJ,Q)S{p-^), 



(102) 



which can be written as a WTI for the vertices with ex- 
ternal (dressed) fermion lines attached: 

Q'^'r^iP, s. A) = x'iP^^Q) - x'iPj^Q) , (103) 

where 



F'^(p, E, A) := S{p + f ) F'"'(p, E, A) Sip - f ) 
x'-ip, Q) S{p + f ) F^(p, Q) S{p - f ) . 



(104) 



We see that the WTI for the IPI part F'''' has the same 
form as the one for the full vertex x^'^ with external 
fermion legs. The remaining equation for Q'^xf''^ is again 
obtained from photon-crossing symmetry. 

Note that the right-hand side of Eq. (101) vanishes in 
the case of the nucleon's Compton amplitude J^": the in- 
verse nucleon propagators for the incoming and outgoing 
nucleon are proportional to the negative-energy projec- 
tors AL, A:^^ and vanish on the nucleon's mass shell, i.e., 
when sandwiched between the positive-energy projectors 
of Eq. (64). This is just the statement that the Compton 
amplitude J^^" — J^" + J^" is transverse: 



Jl-IL' 



0, 



(105) 
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where Jq'^ is the analogue of the Born part (99) for the 
nucleon. The two terms in Eq. (105) vanish individually 
in the Born limit where both intermediate nucleons in 
the Born term are onshell: (P ± E)^ = -M^. In that 
case both nucleon-photon vertices reduce to the onshell 
currents, Eq. (95), and become transverse by themselves. 
However, Q'^ J^'^ also vanishes in general kinematics as 
long as the nucleon-photon vertices in the Born part have 
the specific onshell form of Eq. (95), without positive- 
energy projectors for the intermediate nucleon. The Born 
term constructed in that way is fixed by the nucleon elec- 
tromagnetic form factors and usually subtracted from the 
Compton amplitude. Since the corresponding IPI part 
must be transverse as well, one has a gauge-invariant 
decomposition into 'Born' and 'residual' parts. We will 
return to this point in Sec. VIA. 



C. WTI-preserving Compton vertex 

In the following we want to construct the analogue 
of Eq. (72) for the two-photon case, i.e., the most gen- 
eral fermion Compton vertex that is compatible with its 
Ward-Takahashi identity (97) and free of kinematic sin- 
gularities. Constructions for similar seagull vertices with 
external photon and meson/diquark legs have been de- 
vised in Refs. [58-60]; however, the requirement of Bose 
symmetry for the fermion Compton vertex does not per- 
mit a direct generalization of these studies to the present 
case. The structure of the scalar Compton vertex has also 
been analyzed in Ref. [61], although without a discussion 
of kinematic singularities. 

In analogy to the Ball-Chiu construction, the WTI de- 
termines the fermion Compton vertex up to transverse 
parts with respect to both photons. From inserting (74) 
in (97) we see that it must have the schematic structure 



BC 



1 rprp , 



(106) 



where each term is Bose-symmetric on its own. The first 
two pieces are necessary to satisfy the WTI, whereas Tl^^ 
is transverse with respect to both photons: 

Q'^ r^!^ = , r^!^ Q'^ = . (107) 

The Ball-Chiu part Fg^ is obtained from (72) via 

and fixed by the dressed fermion propagator: 

g'^ F^c = s-\p + S) + s-\p - E) 

The remaining contribution F^^"^ is generated by the 
transverse part of the fermion-photon vertex. It is not 
fully transverse with respect to both photon legs but 
merely satisfies 



In the following we will derive the Ball-Chiu part Fg^ and 
return to the remaining pieces Ti^ and F!^!^ in Sees. VD 
and VE, respectively. 

We exemplify the derivation for a scalar vertex 

T^'{k,Q) = 2k^AB, (111) 

where we have set A{k^) — 0; the generalization to the 
vertex (72) for a spin-Y2 fermion will be straightforward. 
Evaluating the WTIs in that case yields 



Q'^ F^^c = [Ab(p+, Q) + Ab(p7, Q) 
+2p- [AB(p+,g)-AB(p7,Q) 



BC 



-2]9'' 



AB{pt,~Q') + AB{p-,-Q') 
AB{pt,-Q')-AB{pi,~Q') 



(112) 



The difference quotients that appear here were defined 
in Eq. (73). Their differences can have zeros in several 
variables, and our goal is to disentangle those zeros and 
express the remainders in terms of quantities that are free 
of kinematic singularities. The arguments of the function 
B in the last equation are the momenta p ± ^ ± S_ ^ ^nd 
we must study the following momentum dependencies: 

b(p^ + \{Q^ + Q'^)±P-Q±P-Q' ±lQ-Q' 



(113) 

We simplify the notation by abbreviating-'^ 

u^p-Q, u'^p-Q', w^lQ-Q' (114) 

and xa — p^ + \{Q^ + Q'^)- Generalizing from the one- 
dimensional case, we define appropriate difference quo- 
tients via the linear terms (with respect to each variable) 
in a Taylor expansion: 



B{xo + u + u' + w) = Bq + u Au - 
+ uu' A,,,,/ + uw A,„„ + u'w A,,/,, 



u' Au' + w Au 
+ uu'wAuu'v, 



(115) 

Including the signs a, /3, 7 G {±1}, we can write 

Bap-f ■— B{xo + au + I3u' + jw) , (116) 

so that these difference quotients are obtained from the 
inverse relations 





Bo 




"+ + + + + + + +" 




B+++ 


(108) 


uAu 
u' Au' 




+ + + H 

+ + + H 




B++- 
B+-+ 




wAu, 


1 


+ - + -H h- 




B+— 




uu' Auu' 


^ 8 


+ + h + 




B-++ 




UwAuw 




+ -H H h 




B-+- 


(109) 


u'w Au'w 
Uu'w Auu'w 




+ hH h 

+ h- + + - 




B—+ 
B 



(110) 



^ We will use these abbreviations only in the present and the fol- 
lowing subsection; u should not be confused with the Mandel- 
stam variable that we defined in Section IIIB. 
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Inserting Eq. (115) for each instance of the function B 
in Eq. (112) yields 



Q'' V^l = 2 Q^F[ + Ap^iuF!, - wu'F^) , 



where Fi, F2 and F3 are defined by 



UU W 7T 



(117) 



Up to transverse terms in Q'^ and Q'^ , this equation has 
the solution 

V^'c = Gi J^'' + G2 p^'p" + G3 (p^Q'" - Q'Y) ■ (124) 

This is the resulting vertex for a scalar particle; it is 
Bose-symmetric and free of kinematic singularities. If the 
vertex is taken onshell, one has u = u' (because p - A = 0, 
cf. Eq. (22)) and hence Fi = Fi ~ F[, and the d simplify 
to 



F. = 



F. 



A„, - A 



(118) 



and their Bose conjugates {u —u') are given by 



F' - 



u' A„' — w^Ai 
u'^ — 

u' — up' 
A^ 



, A„„. - A 

UU w 



(119) 



We have achieved an intermediate goal: the difference 
quotients defined by Eq. (115) are regular in the limits 
M — >■ 0, u' — > and w — >■ 0; the Fi and F^ are also regular 
in the limits u'^ — !■ and w'^ — )■ w'^; and their prefactors 
in (117) are linear in the respective photon momentum 
Q or Q'. 

To proceed, it will be convenient to work with Bose- 
symmetric combinations of Fi, F^: 



F, := 



F' 



SF, 



u^ — u'^ 



(120) 



The 5Fi are also regular in the limit — w'^. They 
are not all independent: contracting the first equation 
in (117) with Q" , the second with Q"^ and equating both 
yields the relation 

- u^ + u'^ 5Fo 
SFi = F3 + t SFs -uu' — . (121) 
2 w 

Note that from Eqs. (118)-(119), SF2/W is also regular 
in the limit w — >■ 0. If we now define 



Gi = 2Fi + {u- u' f SFi 



Gt = 4F9 - Aw 



F, 



{u - u')^ 



SF. 



w 



SF, 



, (122) 



Ga = 2{u - u') 6F1 , 
we can express Eq. (117) in terms of Gi, G2 and G3 
Q"" r^c = Q"^iGi + u'Gs) + p''{u'G2 - 2w G3) , 



BC 



Qi'iGi - UG3) +p^(uG2 + 2w;G3) . 



(123) 



Gi = 2Fi , G2 = 4 (F2 - w;F3) , G3 = . (125) 

The generalization to the full Ball-Chiu vertex (72), 
which includes also the vector dressing functions Sa and 
A A, proceeds along the same lines. After repeating the 
steps (112) and (115-121) for the fermion dressing func- 
tion A{P), we obtain the final result 

r^c = E i¥ Gf + Gf) rj- + ^ ^Gf rj^ , (126) 



J=5 



where we temporarily define the Bose-(anti-)symmetric 
tensor structures 



flu 



The functions G1...3 are defined in Eq. (122); we have 
attached the superscripts A and B to distinguish the two 
fermion dressing functions they are based upon. The vec- 
tor dressing function A entails four additional functions 
G5...8 which depend on the previous ones: 



G5 = 2Fi , 
-Ge = {u^-u'^)SF,, 
8G7 = {u + u')G2, 
-8G8 = {u-u){G2 + Sw6Fi) . 



(128) 



Altogether there are six independent tensor structures 
associated with G^2 3 and Gf 2 3 which contribute to 



BC 

and (128), 



FgQ. Note that the kinematic prefactors in Eqs. (122) 



u + m' = 2p • E - zZ + yVl - zVl - 
u — u'=p-A^z, 



(129) 



compensate the Bose- and charge-conjugation symme- 
tries of the basis elements, cf. Eq. (57): ti, T2 and T5 
are symmetric (BC = H — h), whereas appropriate prefac- 
tors must be attached to rg (BC — ), tj (BC = — h), 

and T3, Tg (BC = H — ) to ensure invariance of the full 
vertex, Eqs. (54-55). 

To conclude this section, we note that the resulting 
structure of the vertex x^'^ from Eq. (104), where the 
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Born terms are included and external propagator legs 
attached, is completely identical to Eq. (126). The WTI 
for the fermion-photon vertex with external propagator 
legs is given by 



(130) 



which, compared to (70), merely amounts to a replace- 
ment A ^ ay and B — >■ —as, i.e., the dressing functions 
of the inverse fermion propagator are replaced by those 
of the propagator itself: 



1 



A(fc2) fc2 + Af(/c2)2 

M(/c2) 1 



(131) 



The complete analysis carries then through to the Comp- 
ton vertex x'bc- 

D. Partially transverse terms 



We still need to work out the contribution Tip in 
Eq. (106) that is needed to satisfy the Compton WTI 
and depends on the transverse part Tt^ of the fermion- 
photon vertex. With the basis decomposition of Eq. (90) 
one can construct the contributions Ti^'^^ to the Compton 
vertex 



r^''(p,E,A) = ^r^y,(p,E,A) 



(132) 



i=l 



for each basis element in Eq. (89) separately. The right- 
hand side of the Compton WTI will contain differences 
of basis elements as well as sums and differences of 
the dressing functions fi. For the former, we isolate the 
relative momentum fc'^ in each basis element: 

rf(Q) * = 1,3 

Tf(fc,g) = <j k^Tt"iQ) 1 = 4,5 (133) 

fc°fc'3rf"^(Q) 1 = 2,6,7,8. 

For the elements that are linear in k" , the fc— independent 
remainders are 



.M" _ 1 



QQ 



(134) 



where we exploited the definitions in Eqs. (76) and (84). 
For the fc— quadratic basis elements one obtains 



(135) 



where we symmetrized in the indices a and /3. 
The dressing functions 

h{k,Q)^ U{k\k-Q,Q^) 

depend on three Lorentz-invariants, so that the following 
combinations appear in the WTIs: 



h{p^ ,Q) ^ h{p^ + \Q'^ ±u' ,u±w,Q^), 

fi{pf.-Q')^n{p%Q) 

= hip" + \Q''tu,-u' ±w,Q'^). 



(136) 



They are related by Bose exchange Q O — Q', u o —u'. 
In analogy to Eq. (115) we define difference quotients by 

flp ■■ = h {P" + \Q'^ + au',u + I3w, Q2) 

= [fo + au'Au' + I3w Au, + ajSu'w Au'w] ^ , 
(r - /. (P' + IQ' - ^n, -u' + Q'') 

= [/o - auAu^ I3w A„ - aPuw Auw] ^ , 

(137) 

where a, /3 G {±} and /o = /(^'^ + |(^'^ u, Q^). The tilde 
is here merely a reminder that the functions differ from 
those in Eq. (115). The difference quotients are obtained 
from the /^^ via 



(138) 



The following combinations appear in the WTIs: 

h{p^,Q) - MpJ,Q) = /;+ - f-- = Q'^Ht, 

hiP^.Q) + Mpj.q) = /;+ + f-- = m, 

Mpt,-Q') - Mp7,~Q') - if-- - fl+y - Q^H^, 
Mpt-Q') + Mp7,-Q') = if-- + f++y - 2i?,', 

where we can exploit Eq. (137) to obtain 



i7f = (2p^A„, +Q^A^,)., 
iJf = (2p''A„ - Q'"A^)„ 

Hi = ifo + u'wAu'w)i, 

H'i = (/o - uw Auw)i- 







"+ + + +" 




7++" 


v! An' 


1 


+ H 




/+- 


wAyj 


^ 4 


H h - 




/-+ 


u'wAu'w 


i 


.H ^. 







(139) 



Evaluating the WTIs is now straightforward. For the 
two basis elements i = 1,3 that are independent of the 
relative momentum we arrive at: 



Q'^'HtiTl'iQ). 



(140) 
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Since the t!^ are already transverse in their arguments, 
we can simply add both contributions to obtain a Bose- 
symmetric solution of both equations (up to fully trans- 
verse terms): 



-zr^';, = i7rrr(Q) 



For the elements z = 4, 5 we obtain similarly 

+ i/,r7(Q) + ff,'rr(-Q'), 
and for the elements i = 2, 6, 7, 8 we have 



(141) 



(142) 



= iff <(p, Q) + H[\i:{p, -Q') 



- T,i 



\ [i/frr(-Q',Q) +<rf (Q,-g')] (143) 
2A:" [H,r7"(Q) + i7;Tr"(-Q')] • 



Note that for a scalar Compton vertex only Tl^^ con- 
tributes. 

The final result for the WTI-preserving IPI fermion 
Compton vertex is the sum of Eq. (126), which is com- 
pletely fixed by the dressed fermion propagator, plus the 
eight that constitute in Eq. (132). The latter are 
fixed by the transverse part of the fermion-photon ver- 
tex. The full Compton vertex including the Born term is 
then given by 



-r^llf I pM^^ I p/^i^ I pM^ 



(144) 



The first three terms on the r.h.s. are necessary to sat- 
isfy the WTI of Eq. (101). The above construction can 
also be useful in view of the nucleon Compton ampli- 
tude J'"': since the r.h.s. of the WTI vanishes onshell, 
Eq. (144) provides a gauge-invariant decomposition into a 
fully transverse 'generalized Born part' Jg" -f Jg^ -I- Jt^^ 
that depends on the (offshell) nucleon propagator and 
nucleon-photon vertex, and another fully transverse re- 
mainder Jl^'^ which must be determined dynamically. 



E. Fully transverse part of the Compton vertex 

The remaining part F!^!^ of the fermion Compton ver- 
tex in the representation (144) is transverse with respect 
to both Q'^ and Q'^ . Its structure is of particular impor- 
tance for the onshell nucleon Compton amplitude: since 
the r.h.s. of the onshell Compton WTI (105) vanishes, 
both the full amplitude J^'' and the residual part J''" (if 
a suitable gauge-invariant Born term is subtracted) are 
fully transverse and therefore subject to possible kine- 
matic constraints. 

We recall from the analysis in Sec. (IV) that there are 
9 X 8 = 72 independent transverse basis elements in the 
general offshell case and 9 x 2 = 18 transverse basis el- 
ements on the mass shell. In order to study real or vir- 
tual Compton scattering, we must find the analogue to 
Eqs. (89-90). The tensor structures in the resulting basis 



should be free of kinematic singularities in the relevant 
kinematic limits so that the respective dressing functions 
are kinematically independent. For the onshell Compton 
scattering amplitude, such a basis has been constructed 
by Tarrach [62], with a later modification in Refs. [49, 50]. 
In the following we attempt to build such a basis from a 
systematic point of view. 

Construction of the basis. One can generalize the 
construction of Eqs. (84-89) to the Compton case in the 
following way. We start with the ten basis elements that 
are independent of the photon momenta: 



Pi 




aP 
P2 




aP 

Ps 




aP 
Pi 




aP 
P5 




aP 

Pe 




aP 
P7 


^pa^P _^apP^ 


aP 

Ps 




aP 
P9 




aP 

Pw 





7V,3 



(145) 

where the generalized commutator [ • , • , • ] has been de- 
fined in Eq. (76). Second, we want to contract them with 
all possible tensor structures that are transverse with re- 
spect to both photons and free of kinematic singularities. 
To this end, we use the definitions from Eq. (84), 



^ab = 75 

to define the quantities: 



(146) 



plia,Pu, ,N — 1 f na Pi' < t^a Pu 



"'"'{aM--=k{t'Q'a't'bQ±t'Q'b'CQ)^ (147) 



•~fia.p 



where a^,6'^ e {p^, Q^, Q'^}, and primed quantities are 
obtained from 



p, Q, Q' O 



P. Q', Q- 



(148) 



These tensors are symmetric or antisymmetric under 
Bose and charge conjugation, transverse with respect to 
Q'^ and Q"^, linear in all momenta Q, Q', a and 5, and 
free of kinematic singularities. They are not all linearly 
independent; for example, E_(a, 6) and F_(a, &) vanish if 
a — b. 

In the last step we contract Eq. (147) with the p"^ and 
collect the resulting elements in subsets with increasing 
powers in the photon momenta, in order to avoid kine- 
matic relations between the dressing functions. Since 
there are only 72 independent transverse elements, one 
will encounter linear dependencies between the possible 
combinations 



FT'^'iaMpfip), 

QT'f^'Mpfip)- 



(149) 
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The detailed analysis yields a maximum set of 16 inde- 
pendent elements which carry two powers of the photon 
momenta; they can be chosen as 



(150) 




There are 40 further cubic elements 

X pr!6(p), 



± ■ 



{p,Q') 



ip,Q') 



(151) 



and 16 that come with four powers of Q, Q': 



-fia.l3iy 
' + 



{Q,Q)pf..s{p), 
GT'"'{Q\Q)pf.Ap)- 



(152) 



This forms a complete 72-dimensional basis for the fully 
transverse part of the vertex which is free of kinematic 
singularities and constraints. 

On the mass shell, only 72/4 = 18 basis elements re- 
main by virtue of the positive-energy projectors defined 
in Eq. (64). A possible choice for a linearly independent 
onshell basis is 



X < 



^r''''ip,p) 
{p,p) 



A^'^'ip.Q) 



Fr''^"(p,Q) 



'(Q,Q) 



> X 



pfip) 
pf{p) 



X AL 



(153) 



This subset is particularly convenient since, because of 
the identities 



A^ pf (p) a; ^ -e^PP^ pP Ai^^rK 



one only deals with the Dirac structures 



(154) 



A-^ A* 



Ai75AV, A^A;, a/ 757" a;. (155) 



Except for potentially nontrivial relations induced by the 
Levi-Civita symbol in Eq. (154), sandwiching with the 
onshell projectors should therefore not alter the powers 
in the photon momenta. Consequently, one would have 
eight quadratic, eight cubic and two quartic basis ele- 
ments. 

Scalar Compton vertex. Returning to our spe- 
cial case of a scalar Compton vertex, we note that five 
spin-independent combinations are encapsulated both in 



Eqs. (150-152) and Eq. (153), namely those that are pro- 



portional to 



(5"^ and do not depend on G^^"'^": 



'-Q'P^pQ 



Ef(p,p) 



K''iP^P)-t'QVpQ^ 



Q'Q 



^Q'P ^pQ ' 



r± {P, ^) - 2 y-Q'p ^QQ =■= ^Q'Q' ''pQ j ' 
F+ (Q, Q) — tqiQ, tg^Q , 



(156) 



where we abbreviated the contraction of the indices a 
and l3 by 



(157) 



They describe Compton scattering on a scalar particle 
and are identical to the Tf^g given in Eq. (15) in Ref. [63] : 



rppf 



fTip^p)^ 
2F^'^(p,g). 



(158) 



For comparison with Ref. [63], note that on the mass 
shell one has p - Q = p - Q' , and in Euclidean conventions 
all scalar products pick up a minus sign whereas gf^'^ is 
replaced by —6^^'^ . We also swapped the indices p, ^ v. 

Relation with other bases. Next, we would like to 
establish a connection with Tarrach's basis for the on- 
shell nucleon Compton amplitude [62] which has been 
frequently used in the context of dispersion relations; 
see e.g. [2, 49, 64]. We consider the modified version 
by Drechsel et al. who write the onshell scattering am- 
plitude as (Ref. [49], Eq. (A5)): 

= J2 B,(Q\Q'\ Q-Q',p-J:) Ai T^AV , 

^€J (159) 

j = {1...21}\{5,15,16}. 

In comparison to Ref. [49], our scattering amplitude 
J'"' is matrix-valued since we contract with positive- 
energy projectors instead of nucleon spinors {A]^Ui = Ui, 
Uf A;^ = Uf). The T^'^ in the above equation correspond 
to those in Ref. [49]: they are essentially Tarrach's ba- 
sis elements except for the first four, which have been 
replaced by the scalar structures of Ref. [63] ; and the el- 
ements 15 le were swapped by Tig 20,21 which can also 
be found in Tarrach's paper [62]. 

In principle, the T^" will be linear combinations of 
our basis elements in Eq. (153); however, one might ask 
if our construction from above can also be used to gen- 
erate them directly. To that end we extend Eq. (146) to 
accommodate the 7— matrix: 



Z^'' cP-'^'^Pn'^^P 

07 



75 £■ 



a 7^ 



(160) 



20 





^ '■Q'Q 


^11 


= ^+ (7, ) 

= 2iM P'X (7, Q)~2F% {p, Q) + M'uj+ F*; (7, 7) 




= Ff (Q,0) 




= 2iM F^J- (7, Q) - 2 F^" (p, Q) + M^co- F'^'^ (7, 7) 


-'4 


^ 2fl^ {p,Q) 


-^14 


= 1 (2iMGf (7,Q') - [Fr''"^(7,0'),pf ]) 




= |F5;-''-(0,Q)pf 


"-^17 


= Fr(7,7) 


rpf-il^ 

-'8 

^9 


— ^_ (1, 

= 4iA/(2Ff (7,p) + G^-'(7,Q')) 
+2p.EFf (7,7)-8Ff (p,p) 


-'is 

rr-if-LV 

-^19 


= 4G^''(7,0') 

= 2Ft"''"^(Q,Q)pf 

= 4 Fr '"'^ (P, 0) pf + 4iA/ G^;'^ (7, Q) 

+2zAfa;_ (G'!-'(7,Q) - G>t''i'y,Q')) 
= 4 F!;-^" (p, 0) pf + 4iA/ G^J (7, Q) 

+2a./a;+ (Gr(7,Q)-G'l'^(7,Q')). 



TABLE III: Representation of the basis of Refs. [49] with our definitions from Eqs. (145), (147), and (161). We have abbreviated 
a;± = (Q2±Q'2)/(4A^2). 



and define similarly as before 



Fr''''(7,7) 



Q'7 ■ 



2 V "Q'-r "bQ "Q'b' '-'tQ 



(161) 



('y h) — i Tr^" f^'' ± f^" f'''' 
17j f ; •- 2 \^^Q'7 fcQ Q'fc' ^7Q 



where again 6^ e {p'^,Q^,Q'^}- While these combina- 
tions do not yield any new basis elements, they provide a 
relatively compact representation of the basis in Ref. [49] 
which is given in Table III. 

The authors of Ref. [49] note that in virtual Comp- 
ton scattering {Q' = 0) the elements Is^g.ig decouple 
from the cross section when contracted with the photon 
polarization vectors, and Tg^g, Ti2,i3, ?2o,2i become pair- 
wise identical, thus leaving 12 independent elements. Ex- 
pressed in terms of our Lorentz tensors with their inher- 
ent Q'^ dependence, the origin of this behavior is appar- 
ent in Table IIL In our basis of Eq. (153), the same argu- 
ment entails that in virtual Compton scattering F+((3, Q) 
decouples and, upon contraction with the photon polar- 
ization vectors, f^(p,Q) and f^(p,Q) on the one hand 
and G^(p,Q) and G_(p, Q) on the other hand become 
identical, which leaves again 12 independent tensor struc- 
tures. We note that an alternative representation of Ta- 
ble III in the VCS limit in terms of the electromagnetic 
field-strength tensor has been given in Ref. [51]. 

Forward VVCS. Finally, another interesting case is 
the onshell basis (153) in the case of forward VVCS, 
where Q = Q': since A — Q — Q' — 0, the incoming and 
outgoing nucleon momenta are identical, pf = Pt = p, 
and one has A( = A'. = A+. This implies 



A+75 A+ = , 



A+75pf (p)A+ = 0, 



(162) 



so that the tensor structures with G± vanish as they con- 
tain one instance of 75. In general, only five of the 18 
elements remain linearly independent, and we call them 
for the moment 



1 


^El''ip,p)A+, 




fiiy 
2 


= Ff (p,p)A+, 




fiiy 
3 


= E';"'^''(p,p)A+ 


P4 


fiiy 
4 


^Fl^'^''{p,p)A+ 


a/9 A 
P4 


fiiy 
5 


= F'i"'^''(p,g)A_ 


a/9 A 
^P4, A+ 



(163) 



The coefficients of T^'^ 



are related to the structure func- 
tions fx, /l, gxT and gLT in the forward limit. In the 
literature (see, for example, Eq. (114) in Ref. [2]), the 
two spin-independent tensors are typically chosen as Tq' 

and pf^p^, where Tq' = t'qq/Q^ is the usual projector 
defined in Eq. (78) and — Tg^p^ . One can work out 
the relations 



(Ti 



1 



' 2 — 2 ^2 y ' 1 

P^ Q-^ 



(164) 



with our basis elements above, and the spin-dependent 
structures can be written as 



' 3 — 

' 4 



-m2a+£^;(^a+, 



(165) 



where we exploited the definitions (146) and (160). 
T'^^ and (T3 -I- '^4)^^" are proportional to the two spin- 
dependent tensors in Ref. [2] if we define a covariant 'spin 
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vector' via 



yl 1 

2 4 ^ 



(166) 



7t751 



where cr"^ = -| [7", 7*^] and 7^ = T^"" -f". We used the 
properties of the positive-energy projectors: |Z*A+ = A+ 
and A+757t ^+ = ^+757''^+. 

VI. NUCLEON HANDBAG DIAGRAMS 

We have now gathered sufficient knowledge to proceed 
with the calculation of the nucleon Compton scattering 
amplitude. We focus exclusively on the handbag contri- 
bution from Eq. (12) which is depicted in Fig. 10. Within 
the rainbow-ladder truncation, the scattering amplitude 
depends on the DSE solution for the dressed quark prop- 
agator and the solution for the nucleon's bound-state am- 
plitude from its covariant Faddeev equation. For these 
ingredients we refer the reader to the literature [28, 30]. 
We use the quark-gluon interaction of Ref. [39] with the 
input parameters from Ref. [31], and we work with a 
quark mass to„ = md = 3.7 MeV at the renormaliza- 
tion point ^ = 19 GeV so that the resulting pion mass is 
= 138 MeV. 

The remaining component is the quark Compton ver- 
tex that is calculated from its inhomogcncous BSE (18), 
illustrated in Fig. 2. We compute the full vertex includ- 
ing its 128 components in general spacelike kinematics, 
i.e. for 



p'^,t,a>0 and z, y, Z e 



(167) 



with the Lorentz-invariant variables defined in Eq. (28). 
It turns out that the structure of the basis in Table I 
simplifies the analysis enormously; the details of the cal- 
culation are provided in App. B. 

Combining all ingredients yields the handbag contri- 
bution to the nucleon's Compton amplitude J'^''. Its ex- 
plicit Lorentz-Dirac, color-flavor and momentum struc- 
ture is worked out in App. C. The full Compton ampli- 
tude can be spanned onto a transverse basis according 
to Eqs. (60) and (67) or, alternatively, using the regu- 
lar basis elements from Eq. (153). The corresponding 
dressing functions are Lorentz-invariant and can be ex- 
pressed in terms of the four variables {t, X, Y, Z} de- 
fined in Eqs. (22-23), or the variables {t, z^, t, t'} from 
Sec. IIIB, or any other suitable combination thereof. 



A. Born terms and gauge invariance 

The extraction of nucleon polarizabilities is compli- 
cated by the fact that the handbag part does not respect 
electromagnetic gauge invariance on its own. As a con- 
sequence, it is not transverse with respect to the photon 



momenta and the basis elements in Eqs. (60) or (153) 
are not sufficient to fully describe its structure. Unfortu- 
nately the lack of transversality will also affect the trans- 
verse projection of the handbag amplitude as it does not 
satisfy the analyticity constraints that are manifest in 
Eq. (153). 

According to their definition, the polarizabilities are 
related to the Lorentz-invariant dressing functions of the 
residual part J^'^ of the nucleon amplitude. 



(168) 



where the following specific form of the resonant nucleon 
Born part is subtracted from the full amplitude: 



-A^ r^(-Q')^jv(P + s)r5;,(g)A 

-A^ r5^(Q)5Ar(P-E)r^(-Q')A 



(169) 



Here, Sn{p) = {—i^ + M)/{p'^ + A'P) is the nucleon prop- 
agator with nucleon mass M, and 



zr^(Q)=Fi(Q2)^M + 



AM 



(170) 



is the onshell nucleon-photon vertex from Eq. (95) with 
the Dirac and Pauli form factors Fi and F2. If the in- 
termediate nucleon momenta in Eq. (169) are onshell, we 
have 



Q ( ^ p'^-M' , 2MA+(p) 

Sn [p) > p2 + jvf 2 ' ^+ 



1 



(171) 



and the positive-energy projectors ensure that the terms 
in Eq. (169) are transverse for an arbitrary nucleon- 
photon vertex, so that their sum is transverse in the Born 
limit s = u = 0. However, with the ansatz for from 
above, the sum of both terms for Jg"^ in (169) is also 
transverse in general kinematics so that the remainder 
J^'' is gauge invariant as well. This is in general not true 
for other onshell-equivalent formulations of F^, see [57] 
for a detailed discussion. 

In order to obtain J^'^ from a quark-level calculation 
one would need to implement the full set of diagrams 
from Eq. (8): 



Handbag 



Cat's cars 



(172) 



The third term involves the quark six-point function G 
minus the quark Born parts that we shuffled from the 
third to the first term via Eq. (12). The sum of diagrams 
in Eq. (172) satisfies electromagnetic gauge invariance 
by construction and is therefore transverse. The resid- 
ual part is then obtained by subtracting Eq. (169), with 
nucleon electromagnetic form factors calculated consis- 
tently. Thereby Jg"^ cancels the nucleon pole encapsu- 
lated in Jq^ ^ and the transverse remainder encodes the 
polarizabilities. 

While the direct calculation of Jq' might be feasible 
in the long term, it is certainly beyond our present ca- 
pacities. Under the assumption of handbag dominance. 
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FIG. 6: t— channel meson poles in the quark Compton vertex {left panel) and the nucleon Compton scattering amplitude {right 
panel). At the pion pole, the four-point functions separate into the pion's bound-state amplitude/the nucleon-pion current 
(i.e., the pole residue of the nucleon's pseudoscalar current) and the n'yy transition current. 



a sensible alternative would be the construction of a 
'gauge-invariant completion' of the calculated handbag 
diagrams. If the dynamics of J^'^ are dominated by the 
handbag contribution, one can attempt to add a kine- 
matic term 



•^R ~ ""Handbag ^ ""Kincmatii 



(173) 



whose sole purpose is to restore transversality (while 
maintaining analyticity) , so that all non-transverse con- 
tributions in the sum cancel exactly. Such constructions 
are possible and we will explore them in future work. 
For the remaining part we will focus on a feature of the 
Compton amplitude that is unhampered by such prob- 
lems, namely, the t— channel resonance structure. 



B. Meson poles in the t channel 

We recall from Eq. (19) that the quark Compton ver- 
tex can be written as the contraction of the channel qq 
Green function G with the inhomogeneous term Fq"^. If a 
rainbow-ladder truncation is employed, the last quantity 
simplifies to the quark Born diagrams ■ The Green 
function includes all meson-exchange poles at the respec- 
tive location 



t 



''M 



4Af2 4A/2 
where it factorizes into the structure 



G 



A2 



Go. 



(174) 



(175) 



'M 



Here, \1/m defines the meson's homogeneous bound-state 
amplitude and ttim the meson mass, and Gq is the discon- 
nected product of a dressed quark and antiquark propa- 
gator. Via Eq. (19), these poles also appear in the quark 
Compton vertex, cf. Fig. 6, 



as long as the bound-state wave function xm = Gq^m 
has non- vanishing overlap with the Born part Fb, i.e.. 



"'M77 



J XMb,A)F^''(p,S,A) 



^0. 
(177) 

The quantity J^'^^(I],A) defines the two-photon tran- 
sition current for a meson M with mass tom and given 
J^^ quantum numbers. If it is non-zero, the respective 
pole will appear in the quark Compton vertex and con- 
sequently also in the nucleon Compton amplitude. Since 
these contributions are resonant and the onshell transi- 
tion currents are conserved, they are transverse by them- 
selves and we can extract them without concern for ana- 
lyticity constraints. 

It is interesting that the pole structure of the quark 
Compton vertex can already be read off from the basis 
decomposition in Table I. Consider, for example, a scalar 
two-photon current J|^'^^^(S, A). Since the relative mo- 
mentum p is integrated out, we have only the two unit 
vectors and from Eq. (42) for its construction at 
our disposal, and a complete tensor decomposition con- 
sists of the elements 



(178) 



With Eq. (48), they correspond to the elements X^"^, Xg'', 
Xg", Xj'^ and Xg''. Implementing the transversality con- 
straint with respect to the external photon legs leaves 
only two possible transverse Lorentz structures, namely 
Y^^" and Y^^ from Eq. (60), for the onshell current. ^ 
Since the Yj are orthogonal, scalar poles can only appear 
in the dressing functions associated with these two basis 
elements in the quark Compton vertex, combined with 
those Dirac structures in Table I that also contribute to 
the bound-state amplitude of a scalar meson. The same 
observation holds also for the nucleon Compton scatter- 
ing amplitude when expressed in the basis of Eq. (67). 



— 9~ : 



A2 



(176) 



'M 



Expressed in terms of the regular basis from Eq. (156), they are 



linear combinations of tQ/g and tq/Q 
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FIG. 7: Dressing functions and /4,i(t) of the quark 

Compton vertex in the basis (62) at = 0, Z = and a — 1. 
We plot the O"^ Chebyshev moments in z and y. The scalar 
and pion pole locations are shown by dotted vertical lines. 



Similarly, the only possible structure for a pseudoscalar 
two-photon current is Y^" - e^"'"l^ s^d.f' , cf. Eq. (50), so 
that pseudoscalar poles in the Compton vertex can only 
appear in the dressing functions associated with Y^"^. 
The same analysis for the axialvector case, whose onshell 
current is transverse to all three external legs, entails 



that the tensor structures ¥4^^, 



'6 ' 



'7 ' 



Y^" and Y^" 



contribute axialvector-meson poles. 

The dressing functions /i,i(t) and /4,i(i) of the quark 
Compton vertex in the basis decomposition of Eq. (60) 
are shown in Fig. 7. In principle, these functions de- 
pend on the six momentum variables p^, z, y, t, a and 
Z . We performed a Chebyshev expansion in the angular 
variables z and y and show only the zeroth Chebyshev 
moments. Furthermore, we set = Z = and = 1, 
so that the only remaining variable \s t — A^/(4Af^). 
In contrast to the nucleon Compton amplitude, we can 
solve the inhomogeneous BSE for the vertex at timelike 
(i.e., negative) values of t directly since the only needed 
ingredients are the calculated quark propagator and the 
quark-photon vertex in the complex plane. Both dress- 
ing functions /i,i(t) and /4,i(i) diverge, as they should, 
at the respective scalar and pion pole locations which are 
denoted by dotted lines. With our rainbow-ladder input, 
the scalar meson mass is m^^ — 0.67 GeV and the pion 
mass is — 0.14 GeV. 

In order to check our numerical accuracy we can fur- 
ther extract the residues at these poles and compare them 
with known results from the literature. As discussed be- 
fore, the basis elements involving Y^ contain the pion 
t— channel pole via the 77^77 form factor. The respective 
current reads 



>/^;,(S,A) 



;^F^^^(T,r')e'^''"^S"A^ (179) 



where Ft^^^(t^ t') is normalized to unity at r = r' = 0. In 
the chiral limit one has the exact result (77r77 = V^- For 
non-zero current-quark masses, the transition form factor 
F^r-yy in rainbow-ladder has been calculated in Ref. [65]. 
In order to compare with our Compton dressing func- 
tions, we repeat that calculation without giving details. 
The same phase-space considerations from earlier apply 
here, and the accessible spacelike region in r and r' is 
shown in Fig. 4. The result for the transition form factor 
in the symmetric case r = r' is plotted by the thick (red) 
lines in Fig. 8. 

Next, we want to reconstruct the same form factor 
from the quark Compton vertex and the nucleon Comp- 
ton amplitude. Although we calculate the latter only for 
spacelike t, the pion mass is sufficiently small so that the 
extraction of the pole residue is feasible. For — >■ — mj, 
the Compton vertex reduces to 



r^''(p,s,A) 



«'^(P,A)^ 



(180) 



where is the pion's bound-state amplitude 
^^{p, A) = ^ h^\p^,z) Z75 nip, A) WkU Sab , (181) 

i 

with four Dirac structures that can be chosen as S 
{t, p, 4^, 4^]}. We have stated the color and flavor 
factors explicitly (the ak are the Pauli matrices); with 
our color-flavor convention one obtains in the chiral limit: 
h''^\p'^,z) — B{p^)/fT^, where B{p^) is the scalar dress- 
ing function of the inverse quark propagator [66] . Simi- 
larly, the nucleon's Compton amplitude at the pion pole 
becomes 



J^''(P,E,A) 



A^ + mi. 



J^nn{P,A), (182) 



where JttNN is the residue of the nucleon's pseudoscalar 
current at the pion pole: 



(A2)A^*75AV 



2 



(183) 



The rainbow-ladder calculation of the pion-nucleon form 
factor Gt^nn has been performed by us recently [31]. The 
relations (180) and (182) are illustrated in Fig. 6. 

With the bases of Eq. (62) for the quark Compton ver- 
tex and (67) for the nucleon Compton amplitude, the 
tensor structure of the pion transition form factor is pro- 
portional to Y4: 



Yr = 



V2 



2APy/2atil - Z2) 



(184) 



and the pion pole will appear in the respective dressing 
functions f4^a and F4. If we work out Eqs. (180-184) and 
define 



Hit, a, Z) 



v/2ai(l-Z2) V 4AfV 



t + 



(185) 
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FIG. 8: Reconstruction of F-^^j from the quark Compton vertex {left panel) and the nucleon Compton amphtude {right panel) 
via Eq. (186). The thick red line in each panel is the result of the direct calculation for the 7r77 transition current; the thin 
lines are the extrapolated values at various values of t. 



we can extract the form factor at the symmetric point 
T = t' from the relations 

2h\>{p^z) (^gg) 



Fl{t,a,Y,Z) 

G7rAfAr(A2) 



vifhich are valid in the limit t = t'=>Z = 0, a = At — t 
and t — > — / (4M^). We also set = z = y = in the 
quark Compton vertex and y = in the nucleon Comp- 
ton amplitude. The superscript v denotes the flavor pro- 
jection to the isovector part. For the quark Compton 
vertex, flavor is just an overall factor, namely the square 
of the quark charge matrix in Eq. (C4), so that the isovec- 
tor projection is fl = /i/S. For the nucleon Compton 
amplitude, the Dirac-Lorentz and flavor parts arc entan- 
gled, cf. Eq. (C9), and the isovector projection is 



iFl 



Ff 



= Fl - 



(187) 



The extracted transition form factor is shown in Fig. 8 
for various values of t. For t — > — m^/(4M^) « —0.005, 
the curves should converge toward the result from the 
direct onshell calculation of the current. In the case of 
the quark Compton vertex, shown in the left panel, this 
behavior is nicely realized as long is t is not too small. 
For small values of t, the iteration of the vertex BSE 
becomes unstable close to the pion pole. 

The extraction of Fj^-^^ from the scattering amplitude 
in the right panel of Fig. 8 turned out to be more difficult. 
In order to implement the quark Compton vertex dress- 
ing functions /4_i in the nucleon Compton amplitude, we 
are so far restricted to zeroth Chebyshev moments in the 
variables z and y. The reason is that the real relative mo- 
mentum p in the vertex is shifted to p3 = p + P/3 in the 



Compton amplitude, cf. App. C3, where P = (Pi + P/)/2 
is the average nucleon momentum. The vertex depends 
now on r'^'^(p3, E, A), and p^ is complex because is 
negative. Therefore, we must calculate the vertex for 
complex p§ and complex values of the new variables z 
and y that are constructed from p3, S and A. The first 
task is straightforward, but a Chebyshev expansion in 
z and y is impractical because \z\ and \y\ are no longer 
restricted to the interior of a unit circle. If the angular 
dependencies were small, keeping only the zeroth mo- 
ments would be sufflcient, but this is usually only true in 
certain basis decompositions of the vertex. 

For the purpose of generating the right panel of Fig. 8, 
we transformed the eight Dirac basis elements Ta at- 
tached to ¥4"^ to the basis 

1 , ^, A , A] , 9* , tt] , Wu I] , [f, Vt, A] , 

equipped with appropriate prefactors p- A, p-T, and S • A 
to ensure positive charge-conjugation and Bose symme- 
try. I]( is chosen to be transverse to A and p so that only 
the first four elements contribute to the pion pole, as they 
are the only ones that appear in the pion's bound-state 
amplitude. In this basis the angular dependencies close 
to the pion pole are weak and a zeroth Chebyshev approx- 
imation is justified; however, the same approximation in 
a different basis can lead to large inaccuracies. 

We note that the problem could be cured by a moving- 
frame calculation of the quark Compton vertex, where its 
inhomogeneous BSE is solved directly in the kinematics 
that appear in the scattering amplitude. In that case all 
variables are real and no complex continuation is neces- 
sary, however at the price of a dependence on two further 
angular variables. Such a course might turn out indis- 
pensable for the extraction of polarizabilities where the 
full 128-dimensional basis decomposition of the vertex is 



25 



integrated over. We will investigate this option in the fu- 
ture. 



VII. CONCLUSIONS AND OUTLOOK 

We have provided a theoretical analysis as well as first 
numerical results for the nucleon's Compton scattering 
amplitude in the Dyson-Schwinger approach. We have 
previously derived a decomposition of the scattering am- 
plitude at the quark-gluon level that satisfies electromag- 
netic gauge invariance by construction. However, its full 
numerical calculation is complicated by the presence of 
the six-quark Green function which includes the s— and 
It— channel nucleon resonance terms. 

In this paper we have calculated the nonperturbative 
handbag contribution to the scattering amplitude. It is 
microscopically represented by the quark Compton ver- 
tex, which includes the quark Born terms and the quark- 
antiquark Green function in the t channel that com- 
prises all t— channel meson resonances. Based on a well- 
established rainbow-ladder ansatz for the quark-gluon in- 
teraction, we implemented the Dyson-Schwinger solution 
for the dressed quark propagator, the covariant Faddeev- 
equation solution for the nucleon's bound-state ampli- 
tude, and we solved the quark Compton vertex from its 
inhomogeneous Bethe-Salpeter equation. 

For the purpose of our calculations we have constructed 
a 128-dimensional orthogonal tensor basis for the quark 
Compton vertex and studied its symmetry properties. 
The same basis can be used for the nucleon Compton 
scattering amplitude, where the number of independent 
elements reduces to 32 and further down to 18 if the 
transversality conditions are implemented. Moreover, we 
have constructed a fully transverse tensor basis that is 
free of kinematic singularities in the relevant kinematic 
limits. We have established its relation with bases used 
in the literature, and we studied its simplification in the 
limits of virtual and forward Compton scattering. 

We have also gcncralizccl the Ball-Chiu construction 
to the two-photon case and established the general form 
of the fermion (and also scalar) Compton vertex that is 
compatible with electromagnetic gauge invariance. This 
construction can be used for devising models for the scat- 
tering amplitude at the quark level based on the dressed 
quark propagator, the nucleon wave function, and meson 
pole terms in the one- and two-photon vertices. On the 
other hand, it can also be useful for model building on 
the nucleon level, as it provides a gauge-invariant com- 
pletion of the nucleon Born terms in general kinematics 
with an offshell nucleon propagator and nucleon-photon 
vertex. 

We have extracted the 7777 transition form factor at 
the t— channel pion pole, both from the calculated quark 
Compton vertex and the nucleon scattering amplitude. 
We obtained reasonable agreement with the result from 
a direct calculation. We found that the angular depen- 
dencies in the quark Compton vertex are not negligible, 



and in view of obtaining observables in other kinematic 
regions it might be necessary to solve the vertex in a mov- 
ing frame. We note that an analogous analysis for the 
pion electroproduction amplitude, which is possible with 
the methods devised here, can in principle yield input for 
the experimental extraction of the pion's electromagnetic 
form factor from offshell momenta. 

We plan to continue our studies by extracting the ex- 
perimentally measured polarizabilities in virtual and for- 
ward Compton scattering. This task has not been pos- 
sible so far since the handbag contribution is not gauge 
invariant by itself, and therefore not transverse, so that 
the extracted polarizabilities would not satisfy the an- 
alyticity constraints. In the long term, electromagnetic 
gauge invariance can be restored by taking into account 
the full set of diagrams in the decomposition of the scat- 
tering amplitude. As an alternative one can attempt to 
construct ansatze for the scattering amplitude that in- 
clude the dynamics of the calculated handbag diagrams 
while restoring transversality and analyticity. With that 
in mind, we also plan to investigate two-photon contri- 
butions to electromagnetic form factors that are deter- 
mined from the nucleon Compton amplitude at spacelike 
momenta. 
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Appendix A: Euclidean conventions 

We work in Euclidean momentum space with the con- 
ventions 

4 

P-q = '^Pkqk, P^=P-P, f = p-'i. (Al) 

k=l 

A vector p is spacelike if > and timelike if p^ < 0. 
The hermitian 7— matrices 7'^ = (7^)^ satisfy the anti- 
commutation relations {7'',7''} = 25^" , and we define 

^'"' = -^[7^71 , 75 = -7V7V- (A2) 
In the standard representation one has: 

' \iak J ' ' \0 -1 J ' \^10^' 

where dk are the three Pauli matrices. The charge con- 
jugation matrix is given by 

C = 7^72, C'^ = C^ = C-^ = -C , (A3) 



26 



where T denotes a Dirac transpose. Four-momenta are 
conveniently expressed through hyperspherical coordi- 
nates: 



/ \/l — ^1 — y2 sin (j) \ 

2 



V 



(A4) 



and a four-momentum integration reads: 

00 1 1 2ir 



Appendix B: Inhomogeneous BSE for the quark 
Compton vertex 

1. Inhomogeneous BSE 

In this appendix we provide details on solving the in- 
homogeneous Bethe-Salpeter equation (18) for the quark 
Compton vertex. Reattaching indices and momentum 
labels, it has the explicit form 

[r- Tej^^b, S,A) = / /Caa'/3'/3Xjirg.(p',S, A) , 

i' (Bl) 

r^iP, s, A) = S{p+) r^^(p, S, A) S{p^) 

which is illustrated in Fig. 9. S{p±) with p± = p ± ^ 
denotes the dressed quark propagator from Eq. (71), and 
the rainbow-ladder kernel is given by 

4 2 47ra(fc2) . 

= 3 ^2 — p — t;: lYc,^, np'f} , (b2) 

Here, Z2 is the quark renormalization constant, 4/3 is the 
color trace, T^*^ = 5^'^ — k'^k" is the transverse projector 
with respect to the gluon momentum k = p — p' , and 
a(fc^) is the rainbow-ladder quark-gluon interaction from 
Ref. [39]; see also Eqs. (13)-(14) in Ref. [34]. 

The rainbow-ladder BSE respects the WTI for the 
Compton vertex: contracting the BSE with Q"^ and in- 
serting the WTIs (97) and (103) for the Compton vertices 
F'''' and yields the condition 



r''{p+,Q)-T^p-,Q) 

P' 



(B3) 



which is satisfied since the quark-photon vertex is calcu- 
lated from its own BSE, 

V'ip, Q) = Z2 i-i" + j lC^o.'p-p X"(P', Q) , (B4) 
p' 

and the rainbow-ladder kernel only depends on the ex- 
change momentum p — p' . 



2. ExpHcit form of the BSE 

We work in the frame defined by Eq. (26) for the exter- 
nal momenta S and A, Eq. (29) for the external relative 
momentum p, and 



( VT^^z^v/T^^y^ sin0 \ 
Vl - z'^ Vl^^ cos (j) 



(B5) 



for the internal relative momentum. The unit vectors 
that characterize the external momenta are then given by 
Eq. (51). The unit vectors constructed from the internal 
relative momentum p' are defined in analogy to Eqs. (42) 
and (45): d and s remain unchanged, but we have in 
addition 



Pt 



which yields 



(B6) 



(B7) 



Lorentz-invariant relations that combine external and in- 
ternal basis elements, but otherwise involve no further 
momentum, can therefore only depend on the variable 
r ■ r' = V ■ v' = cos <p. Furthermore, in analogy to the set 
^ — {p^i z, y} of Eq. (28), we collect the loop variables 
via 



/ sine/) \ 




/ cos (j) \ 


cos^ 




— sin (f) 1 


I I ] 




^ I ] 









y ■.= [p'\ z'=p'-A, y'^p'^-Er] 



(B8) 



To start with, we express the vertices F^"^ and x!^'^ 
by their basis decomposition (53), where fi^a and fi^a 
denote the Lorentz-invariant dressing functions in either 
case. Second, we exploit the orthogonality relation (59) 
for the basis elements by multiplying the equations with 
rj'^ from the left and taking traces. This yields: 



/,,a(c^. A) = f^jLJ, A) + / g{k^) K,,^ab Ubi^', A) , 
Jp' 



/j,a(w. A) = Gij^abi^}, t) fj,b{l^, A) , 

where 



(B9) 



ab = iTr {fZir, s, d) Y r^;{r', s, d) 7^} , 



Gtj.ab = jTr |r^'^(r, s, d) S{p+) T^'^{r, s, d) S'(p_)| , 



fc2 



(BIO) 



The simplicity of the rainbow-ladder kernel, which de- 
pends only on the gluon momentum k = p—p' and enters 
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K 



FIG. 9: Inhomogeneous BSE for the quark Compton vertex, Eqs. (18) and (Bl). 



through the projector T^*^ = SP"' — k^k"' , entails that the 
kernel matrix Kij^ab only depends on the external and in- 
ternal relative momenta p and p' . The propagator matrix 
Gij.ab does not involve any integration and thus depends 
only on p, but now also on the external momentum A, 
i.e., on the Lorentz-invariant variables uj = {p^, z, y} and 
t. The remaining invariants a and Z enter the equation 
merely as external parameters in the driving term f^^. 

The typical strategy in an iterative solution of Eq. (B9) 
would be to compute and store the kernel and propagator 
matrices for given values of the external parameters t, a 
and Z in advance, i.e., before the actual iteration is per- 
formed. Upon integrating the combination g{k^) Kij ab 
in Eq. (B9) over the angle (f>, which does not appear in 
the dressing functions, we obtain a kernel of the form 



Kij^abiP^^P ,z,z',y,y') 



k: 



ij.ab 



(BU) 



The additional indices in the second step shall indicate 
that we have employed further Chebyshev expansions in 
the angular variables for efficiency, which is justified since 
the angular dependencies are typically small. Unfor- 
tunately, even with economical numerical accuracy, the 
above expression requires several hundred GB of memory 
for storage, and so we must look for better strategies. 



3. Breaking down the kernel 

Let us analyze the kernel in Eq. (BIO) in more detail. 
We can exploit the basis decomposition (52) for the quark 
Compton vertex and write 



if,,,afc=Xrx;.^'^iTr<|f, 



where we have abbreviated 
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Xr(r,s,d), X 



Xf(r',s,d) (B13) 



and similarly for the r's. We have also contracted the 7- 
matrices with the gluon projector, so that 7^ ^ jP — ^/k^ 
is transverse to the gluon. The simple structure of the 
basis elements, where Lorentz and Dirac parts are com- 
pletely factorized, should have consequences for the ker- 
nel as well. X^'^X^'^ — 6ij is diagonal in i, j, and the only 



complication in X^'^X^^'' is the additional dependence 
on r'. Thus, the only scalar product that can enter the 
Lorentz trace is r • r', which suggests near orthogonality 
as well. 

The following Dirac traces can occur in Eq. (B12): 
1 X 1 -> i Tr {f ,7^ 7^} =: Kab , 



75x1^ iTr{f,75 7^r^7^}=: 
1x75-^ iTr{r,7^75r^7^} = 



K 



ab 



Kb 
(r-r')2 ' 



75 X 75 ^ i Tr {r<j75 7t 75 7t} = -^ab , 



(B14) 



where we have pulled out the prefactor in i^^j^ for later 
convenience. 

Combining this with the Lorentz traces, the result for 
the kernel Kij ^b is given in Table IV. It is a 16 x 16 
matrix in the indices where each entry is a further 
8x8 matrix proportional to Kab or K'^^^. We abbreviated 
c = r ■ r' — cos (j). With the basis of Table I, the kernel 
is indeed almost diagonal in i,j (entries not shown are 
zero). We observe that so far we have reduced the 128 x 
128 kernel elements to 5 x 64 independent entries 



Kab, cKab, C^Kab, C Kf, 



(B15) 



that appear in Table IV and have to be integrated over 
separately. It is interesting to note that this structure 
is quite general since it is merely the consequence of our 
choice of basis. The actual ansatz for the Bethe-Salpeter 
kernel only enters via the Dirac traces K and K' to which 
we will turn now. 

The 8x8 matrices Kab and K'^f^ defined in Eq. (B14) 
depend on the momenta fc, r, r', s and d, and therefore 
only on the Lorentz-invariants 



r ■ r' , k ■ r , k ■ r' , k ■ s , k ■ d . 



(B16) 



Since the normalized gluon momentum enters through 
the projector, k can only appear quadratically. The ex- 
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TABLE IV: Full kernel of the Bethe-Salpeter equation for the quark Compton vertex. The kernels K and K' are given in 
Eqs. (B17) and (B20), respectively. 



plicit calculation yields 
/ ai 

Kab = 




(B17) 



where an empty slot again means that the respective en- 
try is zero. Kat is characterized by 12 independent en- 
tries: 

ai = 3 , 

02 = -{r-r') - 2(r • k){k ■ r') , 

as = - 



-1 



2(fc.s)2, 

04 = -1 -2(fc-d)2, 

05 = -a2 - (fli + a3){r ■ r') , 
(oi + ai){r ■ r') , 



(BIS) 



-a2 



and 



bi 


= -2(s • 


k) {k 
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-2(s • 


k) {k 
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^-2{d- 


k) {k 
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b'2 


= -2(d- 


k) {k 


r'), 


b3 


= -2{d- 
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= -bs [r 
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with ar — — (ai + 03 + 04), ag = —(02 + 05 + ae). The 
opposite-parity kernel becomes: 




, (B20) 



with 



Co 

c'o 
ci 

C3 
C4 



1-Ha2r •r'-H2(fc-r')2 
= 1 -I- 02 r • r' -I- 2 (fc • r)^ 
= &i r ■ 
= fe'i r ■ 
= &2 ■ 



r' - 6i , 
r'-bi, 
r'~b',, 
r' - &2 , 

= Co - (ai + 04) - (02 + ae) r • r' , 
= Cq - (ai + 04) - (02 + ag) r • r' . 
= -(63 + ^4 r-r'). 



&'2r. 



(B21) 



Since C3 — C3 = cq — Cg, there are only 8 independent 
entries here, which yields 20 in total for Kab and K'^^. 
In combination with the 5 independent elements from 
Eq. (B15), we arrive at 12 x 3 + 8 x 2 = 52 entries for 
the kernel Kij^ab and therefore have reduced memory re- 
quirements by a factor ^ 300. 

The structure of Table IV owes its simplicity to our 
chosen basis, in particular, to the definite properties with 
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respect to the columns S and T in Table I. Operations 
that can be performed on the Lorentz indices alone do not 
affect the inhomogeneous BSE (Bl), so that each {S,T} 
subspace entails a decoupled equation. The same is true 
for operations that can be performed on the external mo- 
menta s and d that enter the equation. (The remaining 
momenta v and r can appear both outside and inside the 
loop integral, via v' and r' , and therefore do not decou- 
ple.) This leads to 11 decoupled blocks in Table IV and 
11 decoupled inhomogeneous BSEs which can be solved 
independently. Their only difference is the projection of 
the driving term onto each basis element, i.e., the /f^ in 
Eq. (B9). 

Analogous observations hold for simpler BSEs, for ex- 
ample those for the vector and axialvcctor vertices whose 
twelve basis elements decouple into transverse and lon- 
gitudinal BSEs with respect to the external total mo- 
mentum [30, 31]. Another closely related example is the 
tensor vertex where our analysis can be immediately ap- 
plied and yields six decoupled equations. 

We can perform a similar dissection for the propagator 
matrix Gij^ab in Eq. (BIO). Evaluating Eq. (52), with 



P± = P ± t , yields 



Gij^ab — 5ij J Tr <^ Ta 



1 

75 



. S{p+) 



1 

75 



= \ i:v{TaS{±p+)uS{p^)} . 



(B22) 



Since no integration is involved here, all basis elements 
depend only on the external momenta r, s and d, and 
the Lorentz trace is diagonal in This also means 

that only the diagonal components 1x1 and 75 x 75 sur- 
vive. The resulting two 8x8 matrices in a, b are, however, 
more complicated than their kernel counterparts because 
of the momentum dependence of the quark propagators, 
and it is more convenient to perform these traces numer- 
ically. Gij,ab carries an explicit dependence on the exter- 
nal variable t, and a Chebyshev expansion analogous to 
Eq. (Bll) yields the form 



Gij,abip'^,z,y,t) 



G 



mm .nn 



(B23) 



Appendix C: Nucleon handbag diagrams 

In the following we detail the calculation of the hand- 
bag contribution to the nucleon's Compton scattering 
amplitude, illustrated in Fig. 10. In rainbow-ladder trun- 
cation, it is the sum of the impulse-approximation and 
spectator-kernel diagrams. In principle one must add up 
all three permutations a — 1,2,3, where a is the label 
of the quark that couples to the current. However, the 
symmetry properties of the nucleon amplitude relate the 
permuted diagrams among each other so that it is suf- 
ficient to consider only one of the three permutations 
explicitly. In the following we choose a = 3, which cor- 



responds to a coupling to the upper quark line as shown 
in Fig. 10. 



1. Nucleon bound-state amplitude 

Wc first collect the properties of the nucleon's covari- 
ant bound-state amplitude that appears in the diagrams. 
The nucleon amplitude including its full Dirac, flavor and 
color dependence reads 



*(p,(Z,P) 



sabc 



(CI) 



where \^n]ap-y5{p, 9, P) is the spin- momentum amplitude 
and [F„]a6c(i are the two flavor tensors. They transform 
as doublets under the permutation group so that the 
bracket in Eq. (CI) is a permutation-group singlet. The 
Dirac amplitudes ^„ carry three spinor indices a, /3, 7 for 
the quark legs and one spinor index 5 for the nucleon, and 
they are mixed-antisymmetric (5'i) or mixed-symmetric 
(^'2) under exchange of the indices a, /3 and related quark 
momenta. They can be decomposed in a basis of 64 or- 
thonormal Dirac structures X^^^^ which correspond to 
s— , p~ and d— waves in the nucleon's rest frame; see 
Refs. [28, 30, 67] for details: 



64 



(C2) 



The momentum-dependent dressing functions /" are 
the Lorentz-invariant solutions of the covariant Faddeev 
equation. They depend on the five Lorentz invariants 
that can be constructed from the momenta p, q and P 
[P^ = —M^ is fixed). The dressing functions for n = 1, 2 
are not independent but related via permutation-group 
symmetry. 

Similarly, the two isospin-Y2 fiavor tensors F„ carry 
three isospin indices a, b, c for the quarks and one for the 
nucleon. They are mixed-antisymmetric (Fi) or mixed- 
symmetric (F2) with respect to a, b: 



[Fl]abcd = ^ [«0'2]ab<5cd , 
[F2]abcd = *f2]ab Ccd , 



(C3) 



where the cTfe are the Pauli matrices, and they are nor- 
malized to unity: [Fjj,]bad'c[F„]abcd = ^n'nl^d'd- To projcct 

onto proton or neutron flavor states, the nucleon index d 
must be contracted with either of the two isospin vectors 
(1,0) or (0,1), respectively. Finally, the antisymmetric 
tensor eabc in Eq. (CI) denotes the color part which is 
also normalized to unity via the factor l/\/6- 



2. Flavor 

We proceed by combining the ingredients to the Comp- 
ton scattering amplitude. Starting with the flavor traces. 
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the quark Compton vertex comes with a flavor factor 



9n 





1 



9 
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i.e., the squared quark charge matrix. We have also given 
the isoscalar-isovector decomposition with isospin matri- 
ces Us = 1/2 and Uv = o'3/2, where 0-3 = diag(l,— 1) is 
the Pauli matrix. The rainbow-ladder kernel is flavor- 
independent; hence, the flavor trace in both a — 3 dia- 
grams is given by 



p(3) 



it 



d'd 



I bad'( 



Qc'c [F. 



abed ■ 



(C5) 



with d' = d = 1 for the proton and d' = d = 2 for 
the neutron. We keep the doublet indices n' and n for 
the outgoing and incoming nucleon amplitude general so 
we can treat the Dirac and flavor parts in the scattering 
amplitude separately. The above expression vanishes if 
n 7^ n', and for n = n' it yields: 



p(3) _ 



ll 
I f qI+ 2q, 



p(3) _ ^1 9; + 
- 3 I 



ql 



(C6) 



which act upon the doublet indices n',ri (here without 
summation) via [30] 



t(1) p(l) 

n' n n'n 

7(2) p(2) 

n'n n'n 



Inserting this together with(C6) into Eq. (C7) yields the 
result for the total scattering amplitude: 



where J„„ are the Dirac contributions to the a = 3 dia- 
grams. The first line in (C8) is the result for the proton 
and the second corresponds to the neutron. Correspond- 
ingly, the isoscalar-isovector decomposition reads 



22 



J. 



(3) 



3J, 
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with the isospin matrices from Eq. (C4). 



where the matrices are again defined in isospin space. 
The total scattering amplitude is then the sum over the 
three permutations and the final and initial doublet con- 
figurations: 



J^siP,^,^) 



d'd 



a=l ' 



S'S 



-(a) 



d'd 



(C7) 



where the first bracket includes the Dirac-Lorentz and 
color parts to which we will return below. The permuted 
diagrams for a = 1,2 can be inferred from the a = 3 
diagram by applying the two-dimensional matrix repre- 
sentations A4' , Ai" of the permutation group S^, 



M' 



-1 \/3 

-Vs -I 



3. Dirac structure 

The remaining task is to determine the Dirac-Lorentz 
and color structure of the a = 3 diagram. The color trace 
in the impulse-approximation diagram equals 1 while the 
second diagram involving the rainbow-ladder kernel picks 
up a color factor —2/3. We can combine both diagrams 
and write: 



7(3)' 

^Ji'Tl 



S'S 



[i'n']p'a'S'Y{Pf,qf,Pf) 



X s^-MSp'fM 5(p+)r^"^(p3,s,A)5(p3-) 



7'7 



(CIO) 
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The momentum routing is illustrated in Fig. 10. P, E 
and A are the external momenta according to their defi- 
nition in Section III A, with and Pf the incoming and 
outgoing nucleon momenta, pi, qi and Pf,qf are the in- 
coming and outgoing relative momenta; pi, p2, ps and 
P3 = P3 i A/2 are the quark momenta; and p and q are 
the (real) loop momenta. For a symmetric momentum- 
partitioning parameter 1/3 the relative momenta are ex- 
plicitly given by 



Pf ^p + 



A 

3 ' " J 

and the quark momenta by 
p P 



9/ 



Pi = -q 



P2 = q' 



p 

2 



P 
3" 



(CU) 
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and P3 = p + P/i. T^^'^ is the quark Compton vertex 
and the S are the dressed quark propagators. Here we 
also see that it is sufficient to know the nucleon bound- 
state amplitudes for real values of pj and as long as 
the four-vector A is real, which is no longer the case for 
timelike values t <Q. 

The second diagram with the kernel appears in 
Eq. (CIO) via 



[^^n^]al}^s{Pi,qi,Pi) = I Kaa'fSp'ik) 



X Sa'a" {Pl - k) Sf3'f3" {p2 + k) 
X [^n]a"l3"'fS{Pt,qi + k,Pi) . 



(C13) 



(3) 

is one of the three diagrams in the Faddeev 
equation, it is not necessary to compute it anew in the 
form-factor calculation. If we attach dressed quark prop- 
agators to the third line of (CIO) and denote the result 
by we can define a nucleon-quark four-point function 



[^n']l3a5''y'{Pf, Qf^Pf) [^n]al3i&{Pi, qi, Pi) 



and write Eq. (CIO) in a more compact form: 



(C14) 



J. 



(3) 



S'S 



5(p+)r'^-(P3,s,A) 
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(C15) 

Note that we could have removed the appearance of 
S{pt) Ets well by absorbing the quark leg into the nu- 
cleon amplitude . The quantity W is 'universal' in the 
sense that it only depends on the nucleon's bound-state 
amplitude, and the quark propagator and quark-gluon 
interaction that enter its calculation in the first place. If 
the quark Compton vertex is replaced by a vector, pseu- 
doscalar or axialvector qq vertex, Eq. (C15) yields the 
nucleon's electromagnetic, pseudoscalar or axial current, 
cf. Appendix D in Ref. [30] andApp. Bin [31]. Neverthe- 
less, we recall that Compton scattering requires further 
diagrams beyond the handbag structure for a consistent 
and electromagnetically gauge-invariant description. 



4. Calculation of W 

While a brute-force calculation of the integral in 
Eq. (C14) is relatively straightforward, it represents also 
the bottleneck in the numerical determination of the scat- 
tering amplitude. We can break up the integral into 
smaller pieces by choosing a suitable orthogonal basis de- 
composition for the incoming and outgoing nucleon am- 
plitudes. Let us write Eq. (C2) in the form 



["^nUp-isip.q.P) 



E 

kluX 



(CI6) 



with vectors V , R and S defined similarly as in Eq. (42): 

5^ = p?^, R^' = g/, V^' = ^t^o^M pi (^17) 

Orthogonal projection is here understood with respect to 
the nucleon momentum P, so that in the nucleon's rest 
frame the vectors V , R, S and P would reduce to the 
Euclidean unit vectors in analogy to Eq. (51). The co- 
efficients fj^i^^ depend on the five Lorentz-invariant mo- 
mentum variables that can be formed from the three mo- 
menta p, q and P, with P^ — — Af ^ fixed. 

The simplest 64-dimensional orthonormal basis for the 
nucleon's bound-state amplitude is defined by the ele- 
ments 



XklujX 



[R, S, P) = [ta Ffc A^75C]„^ [n r, A+]^, , 

(ClSj 



A:,; 1...4, 



Atj = (1 -I- w /*) /2 are the positive- and negative-energy 
projectors that depend on the normalized total momen- 
tum; T_|_ = 1 and r_ =75; and the four Dirac structures 
are given by 

Fo-1, ri-75f, T2 = fi, T3 = $. (C19) 
The basis elements satisfy the orthonormality relation 

(C20) 

where the charge-conjugated elements are given by 



1 yyk' I' uj' \' \/kluj\ c XX X 

4 ^PaS-y ^a^-yS ~ "kk' Oil' Oww' f AA' , 



X^a/7 = [^"^75 A. F, r,] [A+F; r,] , (C21) 

with Fi = — Fi and all other Fj = F^. The charge- 
conjugated nucleon amplitude is 



[^nWs^ip.q.P) = fklu^X^. 



kluX 



(C22) 



kluX 



The decomposition for the amplitude $„ in Eq. (C14) 
that includes the spectator kernel is analogous. 

When implemented in the form-factor diagrams, the 
momentum dependence on {p, q, P} becomes a depen- 
dence on the incoming and outgoing momenta {pi, qt, Pi} 
or {pf,qf,Pf}. The orthogonal unit vectors must be 
adapted accordingly so that one has {Vi, Ri, Si} and 
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{Vf,Rf,Sf}. Inserting (C16) and (C18) into (C14), the 
function W becomes 



WW 



7(5 



(C23) 



where the fyiii^iy denote the dressing functions of 
and the g]^i^x those of $„, and the kernel is given by 



Kwkk'u:.' = Tr {A;^, % ry n ^ A^} 
Since ta = 1 or 75, it will have the form 



(C24) 



Kxxkk'..' = Sw Ki+l^, + /4fe'L' , (C25) 

where the positive and negative-parity contributions can 
be calculated explicitly: 
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(C27) 



with 



Co -2((l + t)5, 



t6^ 



') 



Ca = Uj' ai ■ Pf , 

Cb^ujbf Pi, 

Cab = Co (flj • bf) -UJUj' {Ui ■ Pf){bf ■ Pi) 



(C28) 



Cab 



-uJ e^"""^^ of b)P^P'^ 



and a,6 G {l^,i?,5'}. 

Going further and exploiting Eq. (C25), we write the 
second line of Eq. (C23) as 



x-± ^ V k''^'^ f"' 

n'n,lV,X / y kk' uiui' J k' I' ui' 



±X dkluiX 



(C29) 
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The dependence in the first line of Eq. (C23) appears 
only in the Dirac structures Tj and Tj, . We can separate 
the Dirac from the Lorentz parts by shuffling the mo- 
mentum dependencies from the first into the second line. 
If we perform the sum over W and the g— integration, 
the resulting Lorentz structures depend on the three mo- 
menta p, P and A and have positive or negative parity. 
They carry zero, one or two Lorentz indices, and as a 



consequence of Eq. (C17) they are transverse to the in- 
coming or outgoing nucleon momenta: 



At: 



and so on, with a, 6 € {V,R,S}. We can now apply our 
findings from the Compton-vertex analysis in Sec. IV A to 
construct orthonormal tensor bases for these quantities. 
In the former case, we had to find transverse basis ele- 
ments with respect to the incoming and outgoing photon 
momenta; here we consider transversality with respect to 
incoming and outgoing nucleon momenta. Let us define 
another normalized set of four-vectors, constructed from 
p, P and A: 



tAf' + iVtP'', 
= \/TTtA^-i\/iP^, 
with 



Pt 



(C30) 



We will use the symbols v and r in this context only 
in the remainder of this subsection; they should not be 
confused with the quantities defined in Section IV A. The 
unit vectors {v,r,di} are by construction transverse to 
Pi = P— and {u, r, df} are transverse to P/ = P+y- 
(Since Pi and Pf are linear combinations of P and A, v 
and r are transverse to both.) Thus, one can immediately 
write down the complete set of basis elements for these 
integrals: 



A 

A^a 

Al 

^ab 



1 



(C31) 



with coefficients that depend on the Lorentz-invariants 
p'^, p - P, p - A and A^. Eq. (C31) applies if the resulting 
tensor has positive parity; for negative parity (induced 
by an odd power of = or 5/ = V/) one has instead 
the general form 



A ^ 

A^a ^ 

Al 







(C32) 



Since, with the exception of dy • d^ 7^ 0, the vectors in 
Eq. (C30) are normalized and orthogonal, the coefficients 
of these basis decompositions are obtained simply by con- 
tracting A'(^ and with each basis element. 

If we work out these relations we arrive at the final 
result for Eq. (C14): 

y = E ^kk\ WM. rx] ,,,, [rx ^k K] ,s ' (C33) 

kk'X 
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where the new Dirac structures are given by 



"o 


= 1 , 


"0 


= 1, 


nf 


= ^ 75 , 


n\ 


= 75^ 






ill 








ni 





(C34) 



Note that one obtains 4x4x2 — 32 independent basis ele- 
ments for the nuclcon-quark four-point function Wgigt^^,, 
as it should be for a spinor four-point function with two 
onshell legs. The coefficients are obtained from the inte- 
gral 



-kk'X 



STT" rf -n'P(k)V(k')V(l)V(l') 
Z^'-'kl'^k'l' ■^n'n,ll',V(k)V(l)\ ' 
IV 



(C35) 



where Cq; = Cqj = S^i, and for fc, Z = 1, 2, 3 one has: 



Cli = (ufc)j • {ai)i 



CL = {uk)f ■ {ai)f , (C36) 



with {uk)i e {v,r,di}, {uk)f G {v,r,df} and a; G 
{y, i?, S"}. Positive parity for all coefficients leads to the 
signs V which we define as 



Vik) 



+1 
-1 



k = 0,2,3, 
k = l. 



(C37) 



In precomputing the quantity W(p, P, A) we have seen 
that it does not depend on the average external momen- 
tum S in the form-factor diagrams, so that its Lorentz- 
invariant dressing functions c^jlf^ are only functions of 
p^, p - P, p - A and A^. This also means that the Lorentz 
frame defined by Eqs. (26)-(27) is not the most econom- 
ical choice since the dependence on the external variable 
Y becomes redundant. We compute these functions in- 
stead in the following frame: 



(C38) 




(C39) 



When reinserting the result (C33) for W(p, P, A) into 
the scattering amplitude (C15), we keep A and P from 
above and introduce the average photon momentum S 
via 



S 

M 



z 



(C41) 



The Lorentz-invariant definitions in Eqs. (22-23) remain 
of course unaltered by this choice. In that way the scat- 
tering amplitude depends on the external variables cr, Y 
and Z only through the quark Compton vertex, p is now 
the loop momentum and we have an additional depen- 
dence on the angle (f>: 



VT 



sm (; 
cos ( 



y 



I 



(C42) 



Note that no interpolation for the coefficients of W is 
required when changing the frame from (C39) to (C42) 
since the explicit form of all variables p^ , p-P, p- A and A^ 
remains the same. With the methods described in this 
appendix we have reduced the original two- and three- 
loop integrals from Fig. 10, which carry a dependence 
on four external variables, to a successive calculation of 
intermediate one-loop diagrams which depend only on a 
subset of these variables. 



5. Extraction of Compton dressing functions 

In the final step we want to extract the Compton 
form factors, i.e., the Lorentz-invariant dressing func- 
tions of the nucleon Compton amplitude J^'^(P, E,A), 
from Eq. (C8). For simplicity we use the basis decom- 
position of Eq. (67) with the dressing functions Fi and 
Gi; dressing functions in other bases can be obtained by 
appropriate basis transformations in the end. If we work 
in the frame defined by Eqs. (C38) and (C41) instead of 
Eqs. (26)-(27), the orthogonal unit vectors v, r, s and d 
for the external momenta assume the form 



with the integration momentum 



q = 



y' sm (f 
y'^ cos( 



I 



(C40) 



The vectors and from Eq. (C30) are then the Eu- 
clidean unit vectors in the 1- and 2-directions whereas 
Vi, Ri, Si and V/, Rf, Sf are more complicated. 






f ' 1 








-Y 




Vl-Y^ 


r = 


Vl-Y^ 




Y 




{ ) 




{ / 



and d = (0,0,0,1) as before. We exploit the orthonor- 
mality relations for the basis elements V^" to write 



1 

75 



(F, + G, ^) A: 



+ ' 



(C43) 
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where one has to pick 75 instead of 1 for the elements 
i = 4, 5, 8, 9. This yields 



t 



i + t 

Y 

VTTt 



Y 
VT+t 



for the elements with i — 1,2, 3, 6, 7 and 



G, = 



1 



2(i + t)(i-y2) 



Tr [^75 Ji] 



(C44) 



(C45) 



for the remaining ones. The functions Fi and Gi depend 
on the four variables t, a, Y and Z. 

To conclude, let us summarize the steps to arrive at 
the Lorentz-invariant dressing functions that constitute 
the handbag part of the nucleon's Compton amplitude: 

• Write the Faddeev-equation solution for the nu- 
cleon bound-state amplitude 5'„, and for the com- 
bination <!>„ that absorbs the spectator kernel, in 
the basis (C16) and extract their invariant dressing 
functions fl]^;, and g^i^^. 

• Combine them in the quantity J- defined in 
Eq. (C29), using the kernels (C26)~(C28). 

• Integrate over q to obtain the Lorentz-invariant co- 
efficients of the four-point function W in Eq. (C35). 

• Contract W of Eq. (C33) with the Dirac struc- 
tures Oss' that are needed to extract the Comp- 
ton form factors. From Eqs. (C44)-(C45) we have 
O S {1, ^, 75, ^75}. Contract the result with 
the precalculated quark Compton vertex according 
to (C15). 

• Add up the permutation-group doublet contribu- 
tions according to Eq. (C8) and extract the Comp- 
ton form factors Fi and Gi for proton and neutron. 

• Apply appropriate transformation matrices to ob- 
tain the dressing fimctions in other bases, for ex- 
ample those defined in Eq. (153) or Table III. 



Appendix D: Basis decomposition into transverse 
and longitudinal parts 

In Section V we have decomposed the fermion Comp- 
ton vertex into a part that satisfies its Ward-Takahashi 
identity and a fully transverse part. Another possible de- 
composition is the one into fully transverse, fully longitu- 
dinal and mixed transverse-longitudinal parts. Both de- 
compositions are not congruent since the WTI-preserving 



contribution has components in all three subspaces. The 
second decomposition is useful for the numerical calcula- 
tion of the quark Compton vertex: when implemented in 
the nuclcon Compton amplitude, only its fully transverse 
part can contribute to physics since the remainders must 
be cancelled by analogous pieces in the non-handbag di- 
agrams of the scattering amplitude. 

We use the usual transverse projectors with respect to 
the incoming and outgoing photon momentum: 



rpfj 



rpfJ 



6"" 



(Dl) 



With their help, the three contributions to the Comp- 
ton vertex, which are independently Bose- and charge- 
conjugation invariant, can be written as 















TT ^ 





pQ;3 



T, 



LL 



Q 
Q' 



2 ^LT 



^ TL 



(02 



(D2) 



where F^^^, F^^, F^^ and T^l have been defined in an 



obvious way. Bose symmetry entails 

F^^(S) = -F^J-E) , Fii(S) = Fii(-E). 



(D3) 



We can now evaluate the WTI, either Eq. (97) for the 
IPI vertex or Eq. (101) for the full vertex, to obtain: 



TL 



Fll = Wi' 



(D4) 



where and — — PFf (— S) generically denote the 
right-hand side of the respective WTI. This yields 



^ LL 
pi/ 

LT 
-pfj. 
^ TL 



W2-Q' 



(D5) 



and therefore all three parts are completely determined 
from the fermion-photon vertices. Note that the M^/* van- 
ish for the onshell Compton amplitude, cf. Eq. (105), so 
that there only the fully transverse piece survives. 

The basis elements in Table I can be rearranged to 
correspond to either one of the three terms in Eq. (D2). 
Since the basis factorizcs the Dirac from the Lorentz 
structure, it is sufhcicnt to analyze the transversality re- 
lations at the level of the 16 Lorentz structures. To single 
out the transverse part F^!^, we evaluate the condition 



g/M pM-^ ^ n Q"" F^" 
With the auxiliary variables 



XZ, b = ^/X\/l- Z^ , 



(D6) 



(D7) 
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where X was defined in Eq. (23), and the orthonormal 
momenta from Eq. (42) we can write 



Q 



^ \ = Vt(fes^ + (aTl)d^) 



so that the transversahty conditions become 
ihs" + {a+l)(r)T^''' = Q. 



(D8) 



(D9) 



The resulting relations between the dressing functions 
allow to construct a fully transverse orthonormal Lorentz 
basis which consists of 9 elements and is given in Eq. (60). 
We note that the 9 x 2 = 18 transverse Lorentz-Dirac 
basis elements for the onshell amplitude, cf. Table II, are 
linear combinations of those in Eq. (153) since both cover 
the same subspace. 

The remaining longitudinal Lorentz structures are con- 
structed analogously. One basis element corresponds to 
the fully longitudinal part and is obtained from 



T/if r""^ = n r^^T/?" = o , 



(DIO) 



namely: 



Yio = 



(1 - X)(X5 - Xe) - 2b (Xg + aX7 + fcXg) 



where we used again 

ni = l + X , 



712 = ni 



(Dll) 



4a^ 
ni 



The remaining seven linear combinations satisfy 



T^? T§-^ = n Q'" r"*^ (D12) 



and are given by 



Yii 
Y12 

Yl3 
Yl4 
Yl5 

Y16 



(l-X)X7 + 26(&X7-aX6) 



(1-X)X8 + 26X5 

(Xi5 + 0X13 + 6X9) 



ni 

= {X,, + aX,, + bXn)-^^ 
1 



ni 
1 



(X16 + aXi4 + 6X10) , 



n2 



(Xi4 + aXi6 + 6X12) - 



2a Y 



15 



^/nrn2 



The complete 128-dimensional basis in Eq. (62), con- 
structed from the 16 Y^ and the eight Tq, is again or- 
thonormal in the sense of Eq. (59). 
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